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ABSTRACT 


Local heat transfer and pressure coefficients around 
a right circular cylinder were experimentally determined at 
Reynolds numbers ranging from 57,000 to 495,000. The turbu- 
lence intensity of the free stream was approximately 0.5%. 
The cylinder was heated externally with a constant heat flux 
and simultaneously cooled in a cross flow of air. A uniform 
heat source was provided by aneweiecas Nichrome ribbon with 
a constant, measured current. The steady state surface tem- 
perature increase above the air temperature was indicated by 
thermocouple information. Employing Newton's law of cool- 
ing, heat transfer coefficients were determined as a func- 
tion of the angular location from the forward stagnation 
line. The free stream dynamic pressure and surface pressure 
distribution around the cylinder were obtained using static ~ 
pressure pickups and their associated manometers. A com- 
parison is made between the heat transfer and pressure data 
collected. Experimental results compared within about six 


per cent of the theoretical solution of Schuh. 
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NOMENCLATURE 


English Letter Symbols 


Ai 


Ao 


Cross-sectional area of wind tunnel 


in settling chamber 


Cross-sectional area of wind tunnel 


at entrance to test section 


Cross-sectional. area of Nichrome 
ribbon 


External surface area of an ele- 
ment of Nichrome ribbon 


Pressure coefficient as a function 


of angular location on cylinder 
surface 


Diameter of cylinder 


Root mean square voltage of hot 
wire anemometer 


subscript; Quantity evaluated at 
film temperature 


Froessling number, F = Nu/vRe 


Convective heat transfer 
coefficient 


Local convective heat transfer 
coefficient as a function of 
angular location 


Current 


Thermal conductivity of air 
evaluated at film 
temperature 
Thermal conductivity of Nichrome 
Nusselt number, Nu = hD/k 


Static pressure of air in 
settling chamber 


ft? 
fii 
Pee 


res 


ft 


volts 


BTU 

hr-ft?-°F 
BTU 

hr-ft?-°F 


amps 


BTU 
hr-ft-°F 
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English Letter Symbols 


P2 


P; 


APi-2 


5 cylinder 


et) 


qy 
Wideal 


actual 


rad 


Re 


Statice pressure of air at 
entrance to test section 


Static pressure of air in test 
section 


Surface pressure on cylinder as 
a function of angular location 


Prandtl number 


Static pressure drop between 
settling chamber and entrance 
to test section 


Pressure difference between sur- 
face pressure on cylinder and 
static pressure at entrance to 
test section, Se sender S Pas 


Heat flux 


Local heat flux from an element 
at some angular location 


Dynamic pressure at any location, i 


Ideal dynamic pressure in test 
section 


Actual dynamic pressure in test 
section 


Heat transfer rate 


Heat transfer rate through 
Nichrome element by conduction 


Heat generation rate within any 


Nichrome element, i 


Heat transfer rate by convection 
from external surface of 
Nichrome element 


Heat transfer rate by radiation 
from external surface of 
Nichrome element 


Reynolds number, Re = UD/v 


Tose ee, > 
ali aia 


in of water 


em of water 


em of water 


BTU/hr-ft? 


BTU/hr-ft? 


lip iv/ fata 
lip y/fatia 
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English Letter Symbols 


R 
au 


15 


AT 


Tu 





Resistance of main heating circuit 

Surface temperature of cylinder 

Local surface temperature of 
cylinder at some angular 
location, 8 


Temperature of ambient air in wind 
tunnel 


Surface Temperature at any element 
jh WalClaheroimlS eablolovoist5 al 


Local temperature difference 


between cylinder surface and the 
ambient air, AT = Ty-T. 
Turbulence intensity level 
Free stream velocity 


PReewst ream VeLocety an Settling 
chamber 


Free stream velocity at entrance 
to test section 


Free stream velocity in test 
section : 


Mean free stream velocity 


Mean fluctuation in free stream 
velocity 


Voltage in heating circuits 
Initial probe voltage of hot wire 
anemometer at zero free stream 


velocity 


Arce “engin onan, element, om 
Nichrome ribbon 


Tunnel Correction Factor 


ohms 


eeu 


eu 


ft/sec 


ft/sec 
ft/sec 
ft/sec 


ft/sec 


ft/sec 


volts 


volts 
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Greek Letter Symbols 


€ 


sb 


wb 


Thermal emissivity of Nichrome 


Wind tunnel boundary correction 
for solid blockage 


Wind tunnel boundary correction 
for wake blockage 


Total wind tunnel boundary correc-— 
tion, €, = Seb Cab 
Density of air in settling chamber 


Density of air at entrance to test 
section 


Density of air in test section 


Uncertainty interval of quantity i 


Kinematic viscosity of air evaluated 


at film temperature 
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Local heat transfer and pressure coefficients around a 
cylinder were experimentally determined at ten different 
Reynolds numbers that ranged from 57,000 to 495,000. The 
turbulence intensity level of the free stream was determined 
with a hot wire anemometer and found to be approximately 
0.5 per cent. 

The cylinder was heated externally with a constant heat 
flux, and simultaneously cooled in a crossflow of air by 
forced convection heat transfer. 

A uniform surface heat flux was provided by helically 
winding Nichrome ribbon into shallow surface grooves that 
extended around the circumference of the cylinder. Approx- 
imately half of the thirty-two inch length of the cylinder 
was wrapped with the main heating circuit. Two smaller 
guard heater circuits were similarly installed on each end 
Ompuae main circuit.  Lhermal energy dissipated by Joulean 
heating, as described by Wood in Ref. 1, was then applied to 
the surface by energizing the electrically resistant ribbon 
with a constant, measured current. The guard heater circuits 
were separately controlled. 

The 4.5 inch diameter cylinder spanned the entire height 
of the wind tunnel test section, permitting a two dimension- 
al study. To prevent any internal free convection, the 


eylinder was filled with foam rubber insulation. 
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After the thermal capacity of the cylinder was saturat-— 
ed, a steady state condition was achieved. Then all of the 
heat added was removed by convection, less any losses due 
to thermal radiation to the tunnel walls or circumferential 
conduction through the Nichrome ribbon. The guard heaters 
barred heat from being conducted out the ends. The local 
surface temperature increase above the ambient air tempera-— 
ture, delta T, was indicated by copper-constantan thermo- 
couples welded to the inner surface of the Nichrome ribbon. 
The test section, which was located at the axial midpoint of 
the cylinder, was instrumented with four of these thermo- 
couples welded to the center wrap of the Nichrome ribbon at 
ninety degree intervals. Since the cylinder was mounted on 
a turntable, delta T could be measured at any desired angular 
position. Employing Newton's law of cooling, local heat 
transfer coefficients were then determined as a function of 
angular location from the forward stagnation line. 

Two static pressure taps were also installed on the 
cylinder, ninety degrees apart. These were connected to 
U-tube manometers widen were referenced against the static 
pressure. Local surface pressure variations around the 
circumference of the cylinder could then be measured direct-— 
ly. The free stream dynamic pressure was measured by two 
static pressure pickups installed in the wind tunnel, and 
indicated on the “air speed" manometer. The measurement 


system, described in detail in Section III, provided 
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sufficient information to evaluate the pressure coefficient 
at any desired angular position. 

Local experimentally determined heat transfer coeffi- 
clents were compared with the measured pressure coeffi- 
cients in order to correlate heat transfer trends with the 
various flow regimes which developed around the circumfer- 
ence) ofthe) cylinder ss wihe Hocatlonjiom the Laninar).. trans — 
tional and turbulent boundary layers as well as separation 
and its attendant wake were found, as expected, to produce 
definite trends in the heat transfer data. 

The Froessling number, Nu//Re, was determined at various 
angular locations and was compared with the experimental 
and) theoretical results of other researchers. The low free 
stream turbulence intensity maintained throughout the pres-— 
ent investigation led to experimental results that were 
Significantly lower than those previously reported. Experi- 
mental results were found to agree within about six per 
cent of the theoretical model developed by Seinuls for predict- 
ing local heat transfer coefficients on the forward portion 


of a-constant heat flux cylinder. 
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II. BACKGROUND 


Almost immediately after beginning a study of heat 
transfer from a cylinder, two absolute facts become appar- 
ent. First comes a total realization that the heat trans— 
fer mechanism is closely interrelated with the associated 
fluid mechanics. An adequate treatment of the former, 
whether it be analytical or experimental in nature, requires 
a basic understanding of the movement, pressure, and general 
behaviour of the fluid particles which conduct and convect 
thermal energy. 

Secondly, it is certainly true that a vast amount of 
research has been devoted to each aspect of the thermal- 
hydrodynamic interface around a heated cylinder placed in a 
uniform stream. .A mere listing of all that has been pub- 
lished would be voluminous. Instead the approach has been 
taken to review the better known, standard works and some 
that are unique. The details of the many factors involved 


are discussed within the framework of this review. 


A. FLUID FLOW AROUND A HEATED CYLINDER 
A brief description of this flow model is summarized 
in this section based on discussions contained in Refs. 2-5. 
The potential flow of an inviscid, ideal fluid around 
the surface of a cylinder consists of an acceleration from 


zero velocity at a forward stagnation point, to a maximum 
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velocity at ninety degrees on both sides, and finally a 
deceleration to zero velocity at the rear stagnation point. 
The corresponding surface pressure distribution around the 
cylinder for ideal flow is simply an inversion of the veloc-— 
ity distribution. The pressure maxima are located at the 
stagnation points and the pressure minima are located at 
ninety degrees. 

The flow of a real fluid around a cylinder is best 
analyzed in terms of changes that occur within the boundary 
layer that forms near the surface. This thin film of fluid 
which grows from a minimum thickness near the forward stag- 
NavGuOiweiicCr ls Tae result Ot thes adherencev or a viscous fiuid 
to the surface, where a no-slip condition of zero velocity 
exists. As long as this layer is attached to a heated sur- 
face it is the medium into which thermal energy is first 
conducted. The fluid motion in the boundary layer then 
eonvects this heated fluid away from the surface. The 
velocity profile through the boundary layer varies as a 
function of the angular location from the forward stagnation 
- dlawanas Hences the Hene transfer rate from the cylinder sur- 
face becomes a function of the angular location. As this 
angular position increases, the boundary layer grows thicker, 
the fluid velocity is increasingly retarded and less heat is 
removed from the surface. 

A boundary layer may be classified as either laminar 
or turbulent. Flow within a laminar boundary layer is char- 


acterized by fluid particles moving in a smooth, orderly 
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fashion and in layers. A turbulent boundary layer is char 
acterized by individual fluid "packets" flowing in irregular 
patterns within the bounds of a composite motion in the 
direction of flow. A turbulent boundary layer is conducive 
to a higher heat transfer rate. 

As the kinetic energy of the free stream attenuates 
within the slower moving boundary layer, the pressure near 
the surface becomes less favorable to the established flow 
direction. Finally, when the remaining kinetic energy of 
the fluid near the surface cannot overcome the adverse pres— 
sure rise, the boundary layer separates, as a reverse flow 
is) initiated near the surface. | Kddies of, turbulence result 
wherein packets of fluid "scrub" the heated surface, signif- 
icantly increasing the heat transfer rate. The resultant 
wake behind a cylinder is characterized by low pressure 
relative to the pressure on the forward half of a cylinder. 
At Reynolds numbers ranging from about 60 to 5000 [5], a 
regular pattern of shed vortices persists ia ine wake; it 
is known as the Karman vortex-street, named after the 
scientist Theodore von Karman. At the higher Reynolds num- 
bers pertaining to this study the wake is characterized by 
irregular turbulent flow. 

Form or pressure drag is a direct result of this low 
pressure region in the wake -caused by boundary layer separa- 
tion. For the range of flow velocities considered in the 
present investigation (50,000<Re<500,000) the form drag is 


large in comparison to skin friction drag. At higher 
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velocities transition to turbulence occurs prior to separa- 
tion. The separation point moves farther back around the 
eylinder and form drag decreases. 

The characteristics of the boundary layer, the location 
of the separation point, and the characteristics of the 
wake and resultant pressure distribution are dependent upon 


the Reynolds number, 


Re = @, (1) 


where U is the free stream velocity, D is the diameter of 
the cylinder, and v is the kinematic viscosity of the fluid. 
The location of the separation point also depends upon the 
magnitude of the free stream turbulence and surface 
roughness. 

Instability first occurs at the rear stagnation point 
at a Reynolds number of about ten, and half of the total 
drag is then form drag. 

The Karman vortex-street first occurs at Reynolds num- 
bers of about 100. Form drag is then predominant. 

At Reynolds numbers between 1,000 and 100,000 a laminar 
boundary layer separates from the surface at an angular lo- 
cation of 80-85 deg from the forward stagnation line. 

Skin friction drag is negligible in this range. 

Above a Reynolds number of 100,000 the fluid in the 
boundary layer may contain sufficient kinetic energy to 
remain attached for a greater distance and not separate until 


reaching an angle of 130 deg with the flow direction. This 
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was shown experimentally by Fage and Falkner [Refs. 6 and 7] 
to occur, if the boundary layer had become turbulent. The 
critical point of transition from a laminar to a turbulent 
boundary layer moves farther back from the front of a cylin- 
der as the Reynolds number is increased. It generally lies 
in the vicinity of 80-95 deg. Boundary layer transition 
prior to separation is given as a cause of the rapid de- 
crease of the drag coefficient at a Reynolds number between 
200,000 and 500,000. 

At Reynolds numbers above 50,000 the heat transfer down- 
stream of separation is larger than on the front half. This 
difference continues to increase with increasing Reynolds 
numbers. 

Schlichting [5] reviews some previous experimental work 
as follows: Hiemenz measured pressure distribution around a 
cylinder in 1911, determining a separation angle of 81 deg. 
His boundary layer calculations provided a comparable separa- 
tion angle of 82 deg. Extensive experimental data on the 
pressure dietbipucion were later collected by Flachsbart. 
His data exhibited a strong dependency on Reynolds number. 
For subcritical Reynolds numbers he recorded a pressure 
minimum at an angle of about 70 deg and a fairly constant 
SSeS over the downstream portion. At Reynolds numbers 
above the critical value he reported a pressure minimum at 
90 deg and generally a closer agreement with potential flow 
theory. Thom investigated the laminar boundary layer on a 


circular cylinder in 1928 at a Reynolds number of 28,000. 
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Fage's study [6] of the boundary layer in the region of 
separation from a cylinder included Reynolds numbers from 
LOO, 000) tor 350,000. “rage concluded that a*critical point 
on a cylinder exists where the boundary layer begins a 
transition to a turbulent state, and that separation occurs 
just beyond the critical point. The critical point was 
shown to be the inflection point in the normal pressure 
distribution curve that immediately follows the point of 
maximum negative pressure. The latter point was also shown 
to move farther around the cylinder with higher Reynolds 
numbers. 

In Ref. 8, Roshko combined free-streamline theory with 
Karman's theory of the vortex-street, and developed a solu- 
tion for drag that is dependent on one experimental input, 
the shedding frequency. In the process, he defined a uni- 
versal wake Strouhal number which is a function only of wake 
Reynolds number. 

Roshko's study of flow past a cylinder at very high 
Reynolds numbers is described in Ref. 9. He observed that 
a sécond transition" range occurs as the drag coefficient 
"increases in the range, 10®<Re<3.5x10°, from a value of 
about 0.3 to 0.7, and then levels off at the latter value." 
He also observed that vortex shedding occurs at Re>3.5x10°. 

Achenbach [10] describes measurements of local pressure 
and skin friction around a cylinder over a range of Reynolds 
numbers from 60,000 to 5,000,000. He determined the loca- 


tions of separation points and transition regions 
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(separation bubbles), and analyzed the data in relation to 
the flow regions involved. 

In Ref. 11, Coder organized and compiled previous ex- 
perimental and theoretical determinations of the boundary 
layer separation point. He also compared this information 
with his own data collected using hot films. His summary 
included separation points at all Reynolds numbers from 
"creeping" flow to "transcritical" flow. His measurements 
were made in the "transition" (or "supercritical") flow 


regime. 


B. THE EVALUATION OF HEAT TRANSFER FROM A CYLINDER 

.in the past twenty years it appears that the emphasis 
has shifted from experimental to analytical studies. 
Reference 12 leads one to conclude that about three approx- 
imate solution methods for heat transfer in the forward 
half prior to boundary layer separation have been proposed 
every two years. This trend has evolved a crucial need for 
eccunave experimental data to verify the theories. It is 
this motivation which led the author to pursue and refine 
an experimental technique. 

A study of forced convection heat transfer from a cylin- 
der may employ one of two separate approaches. These are 
related to either a constant wall temperature cylinder or a 
eonstant heat flux cylinder. In the ROETEE model, delta T 
is held constant and a measured heat flux is increased or 


decreased as necessary to maintain a constant temperature 
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wall. In the latter model, the heat flux is held constant 
and the resulting delta T is measured as a function of a 
angular location. 

Both approaches culminate with an application of 


Newton's law of cooling 


dy = ba = 1, (AM) ee 
68 


where for a specified surface Ag is the local heat trans-— 
fer coefficient at angular location 9, @ is the heat 
(ICEYAS ISIE IPEIGS « Ae Ale PIS SMIPKAGS ChieSel, Cid » slineiwSs9 5 I is 
the local heat flux and (AT) , LS EAS MOC SWEPe EOS Cone 
perature increase above the ambient air temperature. 
Generally, Ay is presented in the dimensionless form of 
the local Nusselt number 
h,D 
AO i «) 


where D is the diameter of the cylinder and k is the ther- 
mal conductivity of the fluid. 

Experimental and analytical results are often presented 
in the form of the Froessling number, named after N. 


Froessling, which is defined by 


= 5 } C4) 
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1. Some Previous Experimental Work 
a. Ernst Schmidt and Karl Wenner 

Schmidt and Wenner present in Ref. 13 their 
well known results for heat transfer over the circumference 
of a heated cylinder. A constant surface temperature cyl- 
inder was used in their experiment. Specifically, a hollow 
brass cylinder was EULESS. AVE one hundred degrees Centi- 
grade by passing steam vapor through it. A londitudinal 
strip was insulated from the vapor. In fact, the strip 
was physically removed and replaced with a copper element 
that contained an internal electrical heating element. By 
controlling the electrical power into the heating element 
_and making use of comparative thermocouple information, the 
strip was maintained at the same temperature as the cylin- 
der wall. Hence, the heat flux was varied to maintain a 
constant wall temperature. A pressure pickup was also in- 
stalled on the cylinder wall. By rotating the entire cyl- 
inder "about a graduated disk," both heat transfer and 
pressure data were collected at various angles with the air 
flow. It was stated that their data were generally only 
taken on half of the cylinder, but symmetry was checked 
with no "substantial difference" between both halves. 

Schmidt and Wenner published the following re- 
sults over a range of Reynolds numbers from 5000 to 426,000: 
"At Reynolds numbers up to around 100,000 high maxima of 
the heat transfer occurred in the forward stagnation point 


at 0 deg and on the rear side at 180 deg, while at around 
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80 deg the heat-transfer coefficient on both sides of the 
cylinder behind the forward stagnation point manifested 
distinct minimums. Two other maximums occurred at around 
115 deg behind the forward stagnation point between 170,000 
and 426,000. At 426,000, the heat transfer at the location 
of these maximums was almost twice as great as in the for- 
ward stagnation point, and the rear half of the cylinder 
diffused about 60% of the entire heat.... On the rear of 
the cylinder at » = 180 deg the heat-transfer coefficients 
for Reynolds numbers of from 25,000 to 426,000 are greater 
than in the forward stagnation point." Heat transfer a- 
round:.a cylinder was stated to be a strong function of the 
Reynolds number and the angle from the forward stagnation 
mnie. Additionally, their results tor heat transfer at* the 
forward stagnation point were said to agree with Squire's 
solution, described in Ref. 14. 

Although Schmidt and Wenner did not directly 
merer tov curbulence factor, Giedt,; in Ref. "15, stated that 
there was a minimum of turbulence and that their data were 
taken in a "relatively narrow jet of air with a nearly 
uniform velocity Glalisiaiealiaionealopas MY 

bow Rov) rand 

In Ref. 16, Fand describes his work which 
involved a constant wall temperature cylinder immersed in 
a uniform stream of water. The Reynolds numbers studied 
were between 11,000 and 64,000. His cylinder was composed 


of a copper core with a surface covering of chromium. 
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Thermocouples and electrical heaters were encapsulated in 
hypodermic containers and embedded longitudinally into 
drilled holes in the test cylinder. The high thermal con- 
ductivity of the copper core enabled a uniform surface tem- 
perature to be faatien estar when the heaters were energized. 
A major design limitation was identified as the low power 
capability of the heaters which provided a temperature in- 
crease of only 4 to 10 °F. Guard heater sections similar 
in design to the test cylinder, were controlled separately 
to prevent any heat flow out the ends. 

Fand's results were shown to be in close agree- 
ment with the McAdams correlation, which had been based 
upon Reynolds numbers from 0.1 to 200. Fand then concluded 
that this correlation which is used to calculate the aver- 


age Nusselt number, 
0-52 0.30 
(Nu), = [0.35 + 0.56(Re), J(Pr), (5) 


could be extended to a Reynolds number of 100,000. Fand 
built a strong argument for a second preferred form of cor- 
relation formula originated by Douglas and Churchill in 


Recay lire 


(Nu) - = a(Re). + b(Re) : (6) 


In this form the first term was said to represent the heat 
transfer through the laminar boundary layer on the front 
portion of the cylinder. The second component, according 
to Fand, represented the heat transfer in the separated 


Region am the rear half. 
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Fand mentioned that in Ref. 18, Van der Hegge 
Zijnen proposed a similar form, and that Richardson in 
Ref. 19 stated that the term representing heat transfer in 
the separated region is a "function of (Re)e ." Fand sum- 


marized these developments and arrived at the following 


correlation formula, 
0,50 0.67 2 
(Nu) » = a(Re ) , + b(Re) p 3 (7) 


He stated that Perkins and Leppert [20] correlated their 
data for water and ethylene glycol with formula (7). How- 
ever, Fand correlated his data successfully with the fol-= 


lowing formula, similar in form to (7), 
0.50 0,58 0,30 3 
(Nu), = [0.35 + 0.34(Re), + 0.15(Re), ](Pr), (8) 


where the term 0.15(Re) represents the heat transferred 
to the separated region. Fand opined that this formula was 
more "related to the physical processes involved" than 
McAdams' formula, and did, in fact, represent ile data 
equally well. Fand further stipulated that equation (8) 
was only valid over one range of Reynolds numbers between 
200 and 100,000. 
Creme rile Gahe dit 

Giedt [15] experimented with a constant heat 
flux cylinder. The relatively large cylinder (4 inch 0O.D.) 
that he used permitted the investigation of Reynolds num- 
bers from 70,800 to 219,000. The turbulence factor of the 


wind tunnel used by Giedt was given as 2.25%. The 
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constant heat flux surface boundary was generated by elec-— 
trically heating a thin Nichrome ribbon that was wrapped 
helically around the cylinder. The temperature rise of the 
cylinder surface above the free stream air temperature was 
measured with thermocouples that were welded to the inner 
surface of the Nichrome. The heating circuit was composed 
of five turns of ribbon, the center turn being the test 
section and the two on either side functioning as guard 
heaters. By maintaining constant power, and measuring the 
temperature increase, Giedt determined the heat transfer 


coefficient from equation (2) 


Q 
Be 


; The center of the hollow cylinder was filled 
with glass wool to reduce conduction losses, which were 
determined to be a maximum of two percent. A static-pres- 
sure tap was also installed on the cylinder. Then, by ro- 
tavangs whey cylinder, both heat transfer and pressure 
information were collected at any desired angle to the flow 
direction. It was stated that in general, data were taken 
on one half of the cylinder circumference but that several 
measurements were made around the entire cylinder circum- 
ference with no substantial differences noted. 

Giedt performed a heat balance on a differential 
element of Nichrome ribbon and thereby accounted for heat 
loss by circumferential conduction through the Nichrome 


ribbon. 
27 
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As Schmidt and Wenner [13] had found previously, 
Giedt again noted that heat transfer is a function of 
Reynolds number and "location along the cylinder circumfer- 
ence." The other major conclusions of Giedt are reiterated 


here for subsequent reference, as follows: 


1. Maximum values of h, occur approximately ten degrees 
further back aillong “the cylinder circumference than 
the minimum static-pressure values. 

eer Och minimum h, values and pressures move back along 
the cylinder with increasing values of Re. 

3. For values of Re up to around 100,000, maximum he 
values occur at the 0 degree and 180 degree points. 

4, At values of Re = 140,000 and up, hy at the rear 
stagnation point is greater than at the forward 
stagnation point. 

De AcueVaiueSs or Re from 1405000 to 2195000, distinct 
maximum hy values occur between 110 and 115 degrees. 

Oe hore this same range of Re, two distinct minimum h 


9 
values occur. The lower of the two is between 80 and 
93 degrees on the front half of the cylinder. The 
second is between 140 and 150 degrees on the back 
half. 
tf. Also for this same range of Re, inflection points 


occur in the pressure distribution curves at the same 


locations as the minimum hy values just noted. 


It is noted that conclusion 3 above agrees 


identically with Schmidt and Wenner [13]. 
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Giedt made a significant comparison of his heat 
transfer results with the skin-friction and static-pressure 
distributions given by Fage and Falkner [6,7]. In his dis- 
cussion for the runs with Re = 140,000 and 219 ,000 he 
stated that "the first minimum of the Nusselt number curves 
coincide in location with the first inflection points of 
the corresponding pressure curves. These inflection points, 
according to Fage and Falkner, mark the location of the 
minimum of the skin friction. The fact that the point unit 
heat transfer coefficients drop to a minimum and then 
abruptly rise near these inflection points also substanti- 
ates their conclusion that transition from laminar to tur- 
bulent flow in the boundary layer occurred at this same 
location." Giedt also concluded that "the second inflec- 
tion points which occur between 149 and 150 degrees in the 
pressure-distribution curves for these same values of Re, 
coincide with the locations of the second minimum hy, val- 
ues." Giedt stated that based upon Fage and Falkners!' work 
"these second inflection points would also indicate the 
locations of second minimums or ESS values of the skin 
Beletionn: Giedt finally expressed a parallel between the 


location of the peak h, value in the vicinity of 110 degrees 


) 
to the peak in the skin friction curve that rises after the 
first minimum. 
di i. Ce Pericinsemdaa. sands G: ueppert 
Perkins and Leppert experimented with two dif- 


ferent techniques using a constant heat flux cylinder. The 


29 


ifs ha) PY fl 0 f 
- cara ae {eet | 


dood aldo noarsqmon sts msotiangte 
emsassig-olists bas NOrtoLsT-ive ode woe reais 
metb'ala wr BYed, veMlint bas ant yd abet ‘HOLT, 
of (000,029 brs 000,00E = ot lage: neurons chia 
sevInNd tIadmuc tiaeevlt si¢ Is mumbo Lae ez? ode” sont 
Io 8inlag tolspalint sett? add Atdw notisool nt 









,Bintog nottoslial essed? .aevaus sivaaesg grihnogn 
vit Yo notzsool sit aAnsm .seniist bak-eysl of t bac 
Jinu gniog otis tedy Io2it sdT —Molsote? ntde odd, to 
; ned? Sas mumtata 2 ot go7b e/nstalttebo, wet ; 


~fingjedue oats etiilog acizoelinat seed? xasea selt 















Tug 
“my of qantmal movi molttenass sett notsefenss ated¥a 


onse els 3p bepqinvso« soewal yrabaued aad at wol th. tag 


~oelint Stones sifd” gadt Sebufsaos’ oslé Shake 4 


Ae ie Fa: 
el 


ad3 Ni seotgoabiOEl brie OE meewssd Suspe Hokaw aantog ane 
"9h Yo zaulev S486 Geedi tol eevine noltudi<gukbead 
-[sv g@ Gumtitian byceek of? to. snoligvol, e497, aale obis 
w1O" ‘Baecle has onset noqv beaed party Seisie JberG ". 
247 egeoldnd cells Sluow etntog nobioelinld baewes ; 
alae oi %o seulavy ote: no enumbado 590008 6 slot? 
wid “eewted Loifsasy | besestgqxe Yllail? 2bed0 » “go kse 
Vilaloty sf? ar site gf 4aeq ot? Yo ROLGa 
2ns teJia see savy ov": AoljJorg7 aly wild: al voeq ond - 


© 


cuntate ten 
) Ore .tL ,fotsees? .O.0 Js 

Uteons isoqqel Onn anislded 

eft {9 au DR. J aod :fanoo 4 gitey esvpiatioes doers? 


v3 


eylinders in both instances were small stainless steel 
tubes that were mounted in a test channel, heated uniformly 
by their electrical resistance to an applied, direct cur- 
rent, and cooled in a cross flow of fluid. 

Their first technique was especially intrigu- 
ing as it employed a constant heat flux cylinder to deter- 
mine the average heat transfer coefficient. As described 
in Ref. 20, a single thermocouple was located along the 
center-line of the cylinder, and the measured center-line 
temperature was mathematically shown to equal the mean 
temperature of the inside surface of the tube. This tem- 
perature was used to calculate the temperature drop through 
the tube wall and thus to determine the surface mean temper- 
ature of the external surface of the tube. The average 
heat transfer coefficient was calculated from the latter 
surface mean temperature. The cooling fluid was water and 
ethylene glycol. L 

Perkins and Leppert correlated the data from 
this first approach at Reynolds numbers from 40 to 100,000 
and Prandtl numbers enon 1 to 300 using several methods of 
correlation. They were most satisfied with results ob- 
tained using the methods of Sieder and Tate for pipe flow 
and Viet and Leppert for flow around spheres. The varia- 
tion in viscosity across the boundary layer was taken into 
account by means of a viscosity ratio factor. Other fluid 


properties were evaluated at the bulk temperature. 
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The other formulas used to correlate this data 
were of two forms, the one proposed by Richardson [19], 
equation (7), and the one proposed by Douglas and Churchill 
[17], equation (6). Perkins and Leppert reported that both 
formulas, after being modified to include the Prandtl num- 
bern and the wiseosity ratio, correlated their data well. 

It was stated that equation (6) as modified might be limited 
to Reynolds numbers not greater than 100,000. 

Perkins and Leppert also illustrated that the 
latter formulas after modification correlate the data of 
some others for a Reynolds number greater than 10. Fand 
[16], however, reported that equation (7) in its modified 
form predicts average Nusselt numbers that are 54-76% 
higher than the experimental results that he obtained. 

Fand suggested that the "free stream turbulence, secondary- 
flows and boundary-layer effects associated with the rec-— 
tangular duct in the apparatus of Perkins and Leppert are 
the cause of this discrepancy." 

In ‘the second experimental technique, described 
in Ref. 21, Perkins and Leppert measured local heat trans- 
fer coefficients. The surface temperature measurement 
method and cooling fluid were different. This time a ther- 
mocouple was attached to the inside wall of the tube and 
utilized to measure temperature as a function of angle. 

The cooling fluid, water, had Prandtl numbers from 1 to 7. 
A pressure pickup was also built into the cylinder to pro- 


vide a pressure measurement capability. Local heat 
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transfer coefficients were determined at Reynolds numbers 
from 2000 to 120,000. 

Both references reflect the significant effort 
made by the two experimenters with regard to blockage 
effects. Potential flow theory was utilized extensively 
for this analysis. 

An analysis was carried out [21] for predicting 
the heat transfer that occurs prior to separation using 
"integral methods with the Pohlhausen velocity profile in 
the momentum boundary layer and a third order temperature 
profile in the thermal boundary layer." 

Perkins and Leppert found [21] that "on the 
rear part of the tube, in the separated region, the tempera- 
ture underwent rapid and large fluctuations." An oscillo- 
graph was used to record "not only the temperature-time 
dependence at each angle in the separated region but also 
provided an accurate way to determine the time averaged 
temperature on that part of the ‘suhse 2 ae Lesal heat 
transfer coefficients were calculated using the temperature 
rise above the ambient air temperature, the measured, con- 
stant heat flux, and the application of Newton's law of 
cooling. The viscosity ratio correction was employed to 
develop local Nusselt number correlations. 

) The turbulence level was reported to be about 
2.9 per cent. It was demonstrated, however, that the data 
collected using both techniques "appears more in iene with 


a l per cent turbulence level." 
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e. dJderome R. Jarcy 

Jarcy conducted an experimental study of heat 
transfer from constant heat flux cylinders at relatively 
high temperatures. Reference 22 describes his study of flow 
field, drag and heat transfer characteristics of a mica 
cylinder. The cylinder was heated electrically using 
Nichrome ribbon wound around the cylinder. The average 
surface temperatures were from 250°F to 1100°F above the 
ambient air temperature. Jarcy reported that the flow for 
his experiment was subcritical and that separation was 
laminar.. He determined this from the pressure distribution 
around the cylinder. Velocities involved were between 40 
and 160 feet per second. 

Jarcy's conclusions stated that "When a cylin- 
der is heated, the value of the dynamic and total pressures 
in the wake of the cylinder decreases as compared to a zero 
heating rate." 

freee Dee ehenehard son 

Richardson considered [19] the separated flow 
region in the rear of a bluff body immersed in a transverse 
stream. He was particularly concerned with turbulent flows 
at high Reynolds numbers that do not reattach after separa- 
tion. Richardson concluded that heat transfer at the rear 
stagnation point "is proportional to the two-thirds power 
of Reynolds number." He also pointed out that it is known 


from the photographic record of Fage and observations of 
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Roshko that in the separated flow region variation occurs 
in the direction of the axis’ of the cylinder: 

Richardson described an approximate method for 
calculating the convective heat transfer in the separated 
flow region of a bluff body and applied this method to a 
cylinder. This method accounts for the separate physical 
mechanisms of heat transfer that exist in the laminar 
boundary layer and separated regions. His method represents 
heat transfer in the separated region as being proportional 
to the two-thirds power of Reynolds number. He mentioned 
previous observations "that separated flow develops first 
at the rear stagnation point and then spreads forward round 
the cylinder at a decreasing rate with increase of Reynolds 
number until a position about 80 deg from the forward stag- 
nation point is reached." He reported that the initial 
development apparently occurs at a Reynolds number of 10. 

In Ref. 23 he stated that blockage effects are felt within 
the separated region. Richardson compared his results with 
the experimental data of others and reported good agreement. 
ge. H.H. Sogin and V.S. Subramanian 

Reference 24 describes an experimental study of 
mass transfer from naphthalene cylinders in a cross flow of 
air by Sogin and Subramanian. They credited this technique 
to Froessling's work in 1938 and 1940, which made use of 
the analogy between constant-property incompressible heat 
transfer and low concentration incompressible mass transfer. 


The specific technique used in their investigation was said 
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to be the same as that of Christian who measured local 
rates of mass transfer from hollow napthalene cylinders in 
axisymmetric flow. Their measurements spanned a range of 
Reynolds numbers from 122,000 to 342,000. Their results 
before boundary layer separation were compared with the 
approximate analytical solutions of H. Schuh and H.J. Merk 
as well as other experimental works. 
Ze OMe me AViare loan Mem ne Ore uakcall Work 

Because of the complexity of the flow pattern once 
the laminar boundary layer separates or becomes turbulent, 
analytical solutions have not been found that predict local 
heat transfer over the entire cylinder. They generally are 
limited to approximate methods for calculating heat transfer 
in the front region of the cylinder extending from the for- 
ward stagnation line back to about eighty degrees. In 
Ref. 12, Spalding and Pun have compiled a summary, classi- 
fication, and comparison of fifteen theoretical analyses of 
heat transfer in this region for a constant temperature 
cylinder. Included are the analyses of Froessling, Seban, 
squire, Eckert, Schuh and la ehchidd | Spalddnge and) Pur 
emphasized the need for experimental data collected at a 
low free stream turbulence level, to act as a standard of 
comparison and assist in the verification and evaluation of 
the proposed solutions. 

Seban and Chan present in Ref. 25 various methods 
for predicting heat transfer to laminar and turbulent 


boundary layers. The variable wall temperature case is 
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considered and, in particular, the constant heat flux case. 
Lighthill's method for laminar boundary layers with almost 
constant properties is also included in the report and 
"recommended for general application, although in certain 
instances, the error in the predicted heat transfer may. 
attain twenty per cent." 

Reference 26 contains Schuh's method for calculat- 
ing laminar heat transfer on cylinders having either a con- 
stant wall temperature or a constant heat flux surface. 
Schuh compares the former with the data of Schmidt and 
Wenner [13] and the latter with the data of Giedt [15]. 
Good agreement was shown at Reynolds numbers below 100,000 
with deviations otherwise. Schuh reports that these devia- 
tions were pronounced at and near the forward stagnation 
point, with the experimental values up to 20% higher. 

inv Ret 27, oguire Us) exact. solution for forced con— 
vection heat transfer ate eine forward stagnation point of a 
cylinder is contained as an unpublished inclusion. His 
resulting value of Nusselt number is given as 


8 De |e 3 
Nu = a, -F | ; (10) 





where O, is a function of the Prandtl number (a table of 


values of a, is provided), and B, is. defined by 


u 


Mod 
B=. GEE) 
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‘The velocity outside the boundary layer is given 


as 


ENO 
u — 


a ee (12) 


where x is the distance from the forward stagnation point 


and U is the velocity of the uniform stream. Then, 


u , 
e- aed, (13) 


and taking for example, an average Prandtl number for air 


as 0.71, one obtains the result 
( ‘ 
NUL = We SISi5) lw : (14) 


3. The Effects of Free Stream Turbulence on Heat 
Transfer Rate 
It has long been recognized that free stream turbu- 
M¥ence has a significant effect on the heat transfer from a 
heated bluff body in cross flow. Some of the major efforts 
to resolve and account for these effects have been that of 
jemekestin (Refs. 28, 29, 30, 31), P.F. Maeder [29], R.T. 
Wood [31], R.A. Seban (Refs. 32, 33) and W.H. Giedt 
(her. 34). Kestin [26] describes a dual effect, one a 
broad change in the heat transfer rate resulting from the 
change in flow pattern, and the other a possible local ef- 
fect at a given angle. Some of his other conclusions are: 
a. The magnitude of the heat transfer variation due to 
the presence of free stream turbulence intensity 


seems to be larger at lower Reynolds numbers. 
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b. The turbulence intensity local effect seems to be 
an increase in the heat transfer rate. 

e. The primary local effect occurs in the heat transfer 
"across a laminar boundary layer developed under the 
influence of a favorable pressure gradient" and is 
"of the order of 10 per cent." For a laminar boundary 
layer near the forward stagnation line of a cylinder, 
the local effect is more pronounced and is "of the 
order of 70-80 per cent for an intensity of turbu- 
lence of 2.5-3.0 per cent." 

d. The local effect on the heat transfer rate is non- 
existent across a turbulent boundary layer or for 


separated flow. 


Kestin mentions [30] that all laminar boundary 
layers "carry stochastic oscillations whose amplitude and 
frequency spectrum is related to the turbulence intensity 
imeche free stream," and that if these disturbances do not 
amplify, the Tolimien-Schlichting theory of stability 
classifies the boundary layer as laminar. 

| Giedt [34] also found an increase of local heat 
transfer rates due to free stream turbulence, including a 
maximum 33.3 per cent increase over the calculated value at 
the forward stagnation point, when a rope net was placed 
about two feet upstream of the leading adige of the cylinder. 

Seban [32] summarized the effect of free stream 


turbulence as an over-all heat transfer increase due to 
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effects such as a lower Reynolds number for transition to 
turbulent flow which permits higher heat transfer rates, an 
increase of heat transfer rate to the laminar layer, and an 
alteration of the character of the flow in the separated 
region which may ultimately cause a movement of the separa- 
tion point. He also reported that for flow with a laminar 
boundary layer the local recovery factor doesn't vary due 
to free stream turbulence alterations. 

Regarding the nature of turbulence, Seban measured 
a longitudinal component intensity and spectral distribu- 
tion. Near the stagnation point he included a previous 
suggestion that a Goertler type instability might occur. 

More recently Kestin and Wood [31] found that the 
instability near the stagnation line of a blunt body is 
dipeerent from Goertler's. As a result of this instability, 
Kestin and Wood reported that "the flow field acquires a 
system of evenly distributed, concentrated vortices arrayed 
in the axial direction and wrapped around the upstream por- 
tion of the body; their spacing is inversely proportional 


to the square root of the Reynolds number. The factor mul- 


wv 


tiplying Re” depends on the turbulence intensity, Tu, of 
the free stream and has a nonvanishing value even for Tu=0, 
decreasing with an increase in Tu." Therefore these authors 
econelude that "the flow field in the laminar boundary layer 
outside a Stagnation line is not two dimensional.... Thus 


the effect of free stream turbulence on heat and mass- 


transfer rates from the forward portion of blunt bodies is 
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indirect. First, a system of vortices is formed under all 
conditions, and this enhances the transfer rates. The 
latter are further enhanced when the intensity of turbu- 
lence is increased because the spacing is decreased 
thereby." 
4, Correlation Formulae For Forced Convection Heat 

Transfer 

Although there does not exist a complete analytical 
solution for heat transfer around a cylinder due to the 
complexity of the flow in the region where the boundary 
layer is turbulent or has separated, design requirements 
for this information must be satisfied. The answers are 
generally provided by means of empirically determined cor- 
relation formulae. All correlations are for average Nusselt 
numbers and not local Nusselt numbers. Although the present 
WOLkeens Concerned with the latter, this discussion of cor= 
relations is included for the sake of completeness. The 
eontributions in this branch of the study have been numer-— 


ous and tailored for almost every conceivable combination of 


-. flow conditions. Perhaps the best known, and one of the 


Pereticst was the work of Hilpert [35] in 1933 whose data 
included Reynolds numbers from about 500 to 200,000. Unfor- 
tunately, the turbulence intensity of his wind tunnel is 
unknown. 

Richardson in Ref. 36 concluded that Hilpert's 
measurements "do not correspond well with present-day under- 


standing of the heat transfer distribution around a circular 
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cylinder." As mentioned before, Richardson subscribed to 
a correlation formula with a separate term to represent 
‘ heat transfer in the separated region, that is proportional 
to Reynolds number to the two-thirds power. 

In Ref. 37 McAdams has put together a compilation 
of much of the work in this area. 

After a review of previous work, Van der Hegge 
Zijnen presented [18] a new correlation formula for combined 
free and forced convection heat transfer from cylinders to 


air. 
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The overall objective of this study was to experimental- 
ly determine local, convective heat transfer coefficients 
around a right circular cylinder placed in a cross flow of 
air exhibiting a low turbulence intensity (~0.5%). To ac- 
complish this objective, an experiment was designed that 
embodied several basic theoretical concepts. These concepts 
include the governing equations for the flow of air around 
a two-dimensional cylinder mounted vertically in a closed, 
low-speed wind tunnel. Also involved are the relations 
that govern the heat generation technique employed and the 


modes of heat transfer to the environment. 


A. DESCRIPTION OF THE FLUID MOTION 
Figure 1 is a schematic diagram of the system consid- 


Emea.  Locatwon one is) upstream, location two is the 
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Figure 1. Schematic Diagram of the Air Flow System 
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entrance to the test section, and location three is the 
center of the test section. All pressure pickups that 
‘were mounted either in the wind tunnel or on the cylinder 
were within this system. The applicable equations for 
calculating the velocity in the test section (section 3) 


are (see nomenclature for definitions of terms): 
1. Continuity Equation (clear tunnel, without cylinder) 
o AUn Ok W (215) 


YI 


2. Bernoulli's Equation (clear tunnel, without cylinder) 
i ‘ il 

2, a es Oya = Is og ly (16) 

3. Dynamic Pressure at Any Location, i 


U. (17) 


See ee 


1 


4, Pressure Coefficient at Any Angular Location on the 


Cylinder, 8@ 


Peto P 
oF = ee 2 (18) 
5p U3 


It was assumed that the density of air was constant 
throughout the system at any instant of time. This assump- 
tin implies that the flow was incompressible, and since 
density is a function of temperature, it also implies that 
any temperature rise through the system was uniform. 


Specifically, 


by = On.= Cy i : (19) 
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Combining equations (15), (16) and (19) and noting that the 
Aerolab tunnel had a contraction ratio of A,/A, a Wp db, 
it is found that 


U, 


2 
re) 
CEE) = AP, 0.99 | 5 | : (20) 





Ideal test section behaviour may be described as the 
state existing if the static pressure and velocity ait loca- 
tion three are identical to those at location two. The 
static pressure variation across the test section was neg- 
Ile Cw RON mpUIcOOsec SiO lmauMem CN ciol mei e mm lnere hore itjine 
static pressure at location two, Pi» Was assumed to be iden- 
tical to P,. 


The dynamic pressure at location three is given as 
2 
q, a %oU, ¢. is (21) 


ideality, U3; = U2, and the ideal dynamic pressure at loca-— 


tion three may be written as - 


Is deal 


= *pU, ‘ (22) 


Combining equations C20) mandy (22) 


2 AP_» 


e002 a s9om 9 ideal (23) 





is written as a second result. The numerator can be meas- 
ured directly. 
It is now convenient to manipulate equation (18) into 


a more suitable form. Since the static pressure was 
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assumed to be constant, the numerator may be rearranged as 


follows 


Re feeb bear sone. é (24) 


This latter form can be measured directly and may be desig- 


nated as 


Bab = AP (25) 


0 eylinder 


The denominator of equation (18) is sometimes described as 


2 
3 


(26) 


4 


Gactual — 20U 


The ratio of q actual to q ideal is termed the tunnel cor- 
rection factor, X, and is a function of the particular wind 


tunnel. For the wind tunnel used, 


q 
X = cactual _ O99) 04 C27) 
Wideal 


Now, combining equations (18), (24), (25), (26) and (27) 


the pressure coefficient is found to equal 


AP A 
linder 
C = — Cytinaer s 28 
Pp APR ~2 a) 


The free stream velocity is found by combining equations 


(econ (26)mand (27). 9°The Pinal result is 


/2° RP —2 


u,v —— . : (29) 


The Reynolds number may then be calculated directly, 


U,D 


Re = 0 e (30) 
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Since the above results were derived from equations 
one and two for a clear tunnel, blockage effects due to the 
presence of a bluff body in a constrained channel have not 
yet been included. Blockage corrections to the free stream 


velocity and Reynolds number are computed in Appendix A. 


B. THERMAL ENERGY GENERATION AND TRANSFER 

Thermal energy was generated by energizing thin Nichrome 
ribbon with a direct current, I. By means of the uniform 
electrically resistance Ones MemNEChicOmenrnh. seheCctinle@allalemeripy, 
was transformed uniformly into thermal energy. This tech- 
nique was described by Wood [1] as a Joulean heating rate, 


where 
Ql= eR. (331) 


Once generated, the thermal energy was dispersed by all 
three modes of heat transfer. Convection was the primary 
mode , and forced convection in particular was the mode of 
interest. References 38-40 support the assumption made 
that external free convection could be neglected. 

Convection ee cnat to the cylinder was considered to 
be negligible, as the cylinder nas Plea with foam rubber 
which prevented any major circulation of air. 

It is helpful to consider an intuitive account of the 
external forced convection involved in the cooling of a 
heated cylinder immersed in a cross flow of air. At the 
forward stagnation point the cylinder receives the full 


strength of the wind and, as one would expect, is 
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relatively cooler. Proceeding back along the sides of the 
cylinder one would expect a continued decrease of the cool- 
ing effect of the wind and a’ gradual increase in the wall 
temperature of the cylinder. Once the boundary layer 
separates or becomes turbulent it seems logical that the 
eylinder wall would rapidly begin to be cooled again, by 
the turbulent packets of fluid "scrubbing" the surface and 
transporting thermal energy from the cylinder surface. 

The concept employed for quantitatively evaluating the 
external forced convection was the same as that used by 
Giedt [15]. It consists of applying Newton's law of 


cooling, 


= = 8, 2 a, (te) (32) 
Ss 


to the previously discussed expression, 
Oy= 17h. 


(See nomenclature for definitions of terms.) Then the basic 
expression for the local heat transfer coefficient can be 


written as 


I?R 


Be, © AG STs )) : (33) 


Conduction of thermal energy out the ends of the cylin- 
der was prevented by the installation of guard heaters. 
Since data were collected at steady state conditions, radial 
conduction of heat energy into the interior of the cylinder 


was neglected. Circumferential conduction of thermal 
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energy through the’ Nichrome was considered. This transfer 
should vary with position since the temperature gradient 

is a variable of position. Near the forward stagnation 
point, plus or minus about 20 deg, the temperature gradient 
should be small and circumferential conduction should have 
a minimal effect. However, one would expect that farther 
around the cylinder the trend would eventually be a flow of 
heat into a relatively cool region. A net flow of heat 
into an element in this region augments the heat generated 
within the element, causing the wall temperature to be 
hieher than cic would be votherwase. In the vicinity of the: 
separation point the temperature gradient should be steep 
and GhcwetnecuOmeca reunherent ualconduct ion, would. be 
maximized. A net flow of heat out of an element in this 
region of relatively high temperature would tend to indicate 
a higher heat transfer rate than an ideal case without 
eircumferential conduction. 

Radiation of electromagnetic thermal energy from the 
cylinder to the tunnel walls was taken into account and 
Eeeanoadsasmanlogs. Radiation losses were relatively small 
because of the small temperature differences involved, but 
were not considered to be negligible. 

‘a elemental heat balance has been developed for an 
element of Nichrome ribbon to derive the detailed expres- 
Sion required for calculating the local heat transfer co- 


Si mrucment, kh The resulting expression contains all 


9° 


applicable modes of heat transfer. The components of this 
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Figure 2. Heat Transfer for an Element of Nichrome Ribbon 


expression are depicted in Figure 2. (See nomenclature for 


definitions of terms.) The heat balance is first written 


mee Oy a Gh, - Gy - ay Ax = 0. (34) 
It is noted that 

a, = -k,A Ss : (35) 

wk ee = ie a AS 5 (36) 

ee jae) i (37) 

ee, Tece(ra-l ase: (38) 
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and 


1, (LER (FI. (39) 


Substituting into the heat balance one has 





laud Na toi ei i bai Eat 
kA (Nee? =| hip Ae CE te) acid Cain) O30) 


Bune at 
Rearranging 
eae Vv 4 4 
k — + |— - 


Rewriting in a finite difference formulation 


KO An v2 4 4 
hg = TAS am Dering Ae male pce hae te (42) 
A (Da) 


This formula is the core of the computer program in Appen- 
dix B that was written for data reduction. 

The above expression doesn't account ort cereaynamdel 
heating of the cylinder, and for the continuous rapid heat- 
ing of the ambient air by the wind tunnel boundaries and 
impeller at velocities above 125 feet per second. Because 
the former effect tended to raise the unheated cylinder 
above its initial temperature and the latter caused the 
slower responding cylinder to lag the ambient air tempera- 
ture, the two effects labored in opposition. Their com- 


bined, resultant effect was calibrated out of the data as 


shown in Appendix C. 
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IV. EXPERIMENTAL APPARATUS 


A. LOW SPEED WIND TUNNEL 

The experimental study was carried out in the low speed 
Aerolab wind tunnel, located in the Aeronautics Laboratory 
in Halligan Hall, U.S. Navy Postgraduate School. A schemat- 
ic diagram of the wind tunnel is provided in Figure 3. 

The turbulence intensity level of the wind tunnel was 
measured using a hot wire anemometer and its associated 
electronic equipment. The turbulence level ranged from 
0.67 for the low speed run to 0.43 for the high speed runs. 
The details of these measurements and the equipment used 
are described in Appendix D. 

The wind tunnel is powered by a 100 hp electric motor, 
and has a four speed transmission. All experimental runs 
conducted at a Reynolds number less than 350,000 were in 
third gear and the high speed runs were soneneted in fourth 
gear. The tunnel’ was capable of achieving a maximum speed 
of 200 mph with a clear test section. However, with the 
Be tinder in place, the peak speed was only about 172 mph. 

Air speed information was determined by observing the 
pressure drop between two static pressure manifolds which 
indicated static pressure in the settling chamber (Figure 1, 
location one) and near the entrance to the test section 
(Figure 1, location two). This static pressure difference, 


AP\_2 , was read on the water-filled air speed manometer. 
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Schematic Diagram of the Wind Tunnel 
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The manifold for each location was connected to four static 
pressure pickups installed in the wind tunnel walls at that 
location. Setting the zero and reading the air speed 
manometer required a great deal of consistent, meticulous 


effort to obtain acceptable accuracy. 


B. CYLINDER DESIGN AND CONSTRUCTION 

All data for this project were collected with a 4.47 
in (0o.d.) cylinder constructed from a cast acrylic resin 
tube. A second cylinder made from phenolic tubing has been 
designed and partially constructed, and is expected to be 
an improved model. 

As is apparent in the descriptions which follow, the 
design concept of both cylinders was the same as that of 
Giedt [15] with only some physical alterations. 

ie 6Lhe Cylinder Made From Acrylic Tubing 

The test cylinder was 32 in long and spanned the 
entire height of the wind tunnel test section. Its 4.47 
in (o.d.) diameter and 0.30 in wall thickness were the final 
dimensions after being machined into round from an original 
W625 in (o.d.). 

The choice of a relatively large cylinder diameter 
permitted sufficiently large Reynolds numbers for examining 
all flow regimes in spite of a fairly limited air speed. 

A circuit for generating a constant heat flux on 
the surface of the cylinder was installed by helically 


wrapping Nichrome ribbon around the circumference of the 
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eylinder. The Nichrome ribbon was 0.375 in wide and 0.0031 
in thick, and had a resistance of 0.456 ohms per foot of 
length. This particular size was selected primarily on the 
basis of availability. It also possessed a sufficient total 
resistance to accommodate enough current to drive the sur- 
face temperature well above the ambient air temperature, 
even for high speed flows. 

To facilitate the wrapping process and to promote 
a smooth surface, the cylinder was grooved 2.5 threads per 
in, 0.375 in wide, and 0.005 in deep. This left a separa- 
tion of only about 0.025 in between individual wraps, giv- 
ing a close approximation to a uniform heat flux. The ends 
of the ribbon were led through slots from, and back into, the 
interior of the cylinder where they were connected to power 
leads. The ribbon was bonded to the cylinder with RTV-108, 
Silicone rubber adhesive. 

The main circuit consisted of 39 wraps and had a 
total resistance of 20.527 ohms as measured on a Rosemount 
commutating bridge, model 920A. Variations in the resist— 
ance due to an increase in temperature were considered 
negligible over the temperature range involved. This cir- 
cuit was located approximately in the center of the tube, 
beginning 9.125 in from the base and extending up 15.125 
in of the tube's length. 

Two smaller guard heater circuits were similarly 
installed above and below the main circuit. These consisted 


of three wraps of Nichrome ribbon whose ends were led 
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through separate slots. Each pair of guard leads were 
then connected to a separate pair of power leads. This 
arrangement allowed the guard heater circuits to be con- 
trolled individually. As shown in Figure 4, the guard 
heater circuits were laid down in the same continuous 
groove as the main circuit and adjacent circuit slots were 
separated by only 0.5 in. The smaller circuits each had a 
resistance of about two ohms. 

While the guard heater circuits were the primary 
means of precluding heat conduction losses out the ends of 
the tube, the wall thickness at each end was also shaved to 
0.0625 in to minimize conduction losses. 

Temperature information was measured using premium 
grade copper-constantan bead thermocouples made frem 0.005 
in wire. One thermocouple was shielded and located on a 
screen in the settling chamber, where the velocity was a 
tenth of the test section velocity. It therefore provided 
a true measure of the ambient air temperature. It was 
referenced against aheenccu UMC uRLOMandeackied as a refer— 
ence for the remaining eight thermocouples which were in- 
stalled on the cylinder surface in the following manner. 
They were led through drilled holes in the cylinder wall 
and were spot welded to the inner surface of the ribbon just 
prior to the sealing of the ribbon to the surface. These 
thermocouples directly indicated the temperature difference 
between the cylinder surface and the ambient air. Their 


location on the cylinder is pictured in Figure 4. 
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Figure 4. Test Cylinder 
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The test section ae designated as the center wrap 
of the main circuit. It consisted of four thermocouples 
attached at 90 deg intervals. The remaining four thermo- 
couples were used solely for guard heater control. Two 
were located on the center wrap of the two guard circuits. 
The other two were located near the edges of the main cir- 
cuit and provided the comparative temperature information 
required for the guard heater system to function. 

To avoid any contamination or damage to the test 
cylinder, the actual test thermocouples were not calibrated. 
Instead, an identical thermocouple was welded to a small 
strip of Nichrome. This representative test thermocouple 
was then calibrated in the Rosemount constant temperature 
bath system which had a platinum resistance thermometer 
connected across the Rosemount commutating bridge as a 
standard. An ice bath was used for the reference junction. 
Pimee the results showed close agreement with the published 
table for copper-constantan thermocouples over the applica- 
bee Pe npeeee range, the table was used for data 
reduction. 

An additional check was performed on a similar 
representative test thermocouple to insure that no extrane- 
ous EMF, as described by Davenport and others in Ref. 41, 
was generated due to its attachment to a direct-current 
heater. The test consisted of pulsing the small ‘strip of 
Nichrome with direct current and observing the response of 


the attached thermocouple on a Honeywell chart recorder. 
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Only an exponential rise and fall of thermally induced EMF 
were observed. No step diieneeer, which would be indicative 
of an extraneous electrical EMF, were observed. This con- 
clusion was in agreement with Davenport [41] for a bead 
thermocouple making point contact with an enteraar-lels Wey 
heated substrate. 

Another simple experiment was performed to gain 
insight into the relatively larger thermal expansion of the 
acrylic tube in comparison to the Nichrome ribbon. The 
objective was also to determine an upper limit for surface 
temperature, based upon the endurance of the acrylic to 
resist a combination of continuous high temperatures and 
potentially deforming hoop stress which was applied by the 
cinching action of the ribbon. A third objective of the 
experiment was to determine the ability of the R.T.V. 
‘sealant pon Saas high temperatures. The details of the 
experiment are given in Appendix E. The results are brief- 
ly summarized as follows. Deformation is probable in still 
air at temperatures above 130°F, and the R.T.V. held solidly 
at. temperatures up to EGO 

As the R.T.V. cured after the Nichrome ribbon was 
bonded to the final test cylinder, a contraction of the 
sealant occurred, resulting in various peaks and recesses 
in the ribbon and therefore a somewhat lumpy surface rather 
than the desired smooth one. The effects of these flaws 


are discussed later in the thesis. 
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Two static pressure pickups were installed in the 
cylinder wall 0.75 in. below the lower guard heater and 
at 90 deg intervals. These were made of 0.040 in (0.d.) 
stainless steel tubing, and had an orifice of 0.032 in 
(i.d.). In this way, surface pressure was measured and 
input to one side of a U-tube water-filled manometer. The 
manometer was referenced against the static pressure at 
the entrance to the test section, P,. 

After the thermocouples and pressure taps were in 
place, the cylinder was filled with Silastic, S-5370, RIV 
Foam, a low density resilient silicone rubber. The foam 
rubber was easily poured into the cylinder and quickly 
expanded, uniformly filling every crevice. It provided 
support for the fine thermocouple wires. Internal free 
convection was practically eliminated, and the foam was 
an excellent insulator, exhibiting a thermal conductivity 
of 0.026 BTU/hr-ft-°F. Unfortunately, a serious disadvan- 
tage of the foam was its large thermal capacity which 
considerably increased the delay before reaching steady 
state conditions. ‘Long experimental runs were the out- 
come of this objectionable feature. 

The cylinder was mounted on a turntable base 
that could be rotated from outside the wind tunnel. A 
disk, that was marked every five degrees, rotated with the 
cylinder base. Therefore, temperature difference and 
surface pressure could be determined at any desired angular 


location. 
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Figures 5 and 6 show the cylinder mounted in the 
test section. 

2. A Proposed Phenolic Test Cylinder 

A second cylinder made from 4 in (0.d.) phenolic 
tubing has been designed and partially constructed. It is 
based upon the same experimental concept as the first 
cylinder, but has several different and potentially better 
features. 

The coefficient of thermal expansion is not as 
lance aSwuncmacryilaie Tube.) DUGES Closer TO thareo the 
Nichrome. Therefore, it is possible that a sealant would 
Nome reGurced Over the enuire Hengsth of the ribbon, if 
required at all. This would preclude the development of 
surface irregularities. 

The phenolic tubing can withstand much higher 
temperatures for long periods of time without experiencing 
deformation. 

A space for four additional thermsceneles has been 
provided on one side of the cylinder in the test section 
région, Beene in thermocouples spaced every 30 deg. 
This feature will expedite data collection. 

At the time of this writing, the phenolic tubing 
has been machined and grooved to accept the three circuits. 
The rotating base has been modified to accept the smaller 
diameter tube, and electrical design requirements have been 


completed. 
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Figure 6. Closeup View of Test Cylinder 
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C. ELECTRICAL CIRCULTRY AND EQUIPMENT 

All power leads and thermocouple wires were led down 
the interior of the cylinder and out of the wind tunnel 
through a hole in the cylinder base. Each pair of power 
leads for the three circuits was connected to individual 
Lambda model LK 345A FM 60-volt regulated power supplies. 
A Non-Linear Systems Inc., Model 481 digital voltmeter was 
used to monitor the power settings for any of the three 
circuits to the nearest miliivolt. 

The thermocouple wires were connected to a ten posi- 
tion rotary switch: eines posteaons for the thermocouples 
on the cylinder, one position for the thermocouple in the 
settling chamber and the tenth position was used to include 
the ice junction. The output leads from the switch were 
connected to either a Honeywell chart recorder during the 
transient, warming up period, or to a Leeds and Northrup 
millivolt potentiometer for the accurate data collection 
at steady state conditions. 

Figures 7 and 8 illustrate the arrangement of the 


equipment. 
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V. EXPERIMENTAL PROCEDURE 


In the development of a precise experimental technique 
for the test cylinder constructed for this study and others 
that might follow, there were three distinct phases of ex- 
perimentation, each spanning several months of research. 

Initially, the heat transfer characteristics of the 
cylinder were investigated at a single Rewaolas number of 
53,0002) During (this time a large amount of preliminary. 
data were recorded to Bet eriine if the temperature distrib- 
ution over one side of the cylinder was the same as on the 
opposite side so that data could be collected from one side 
only; to ascertain which side of the cylinder was prefer- 
able and would provide the smoother data, particularly in 
view of the surface irregularities mentioned previously; 
and to resolve a preferable pattern of rotation in conjunc- 
tion with the preferred side. 

The second phase of experimenting consisted of an in- 
vestigation of eight different Reynolds numbers, and was 
peteinally intended to provide the final results. This 
"would, in fact, have been the case had the author been 
more experienced in the art of wind-tunnel testing. How- 
Ever, av various times in the process of collecting this 
data, flaws in the experimental technique were discovered 
which detracted from the credibility of the results. For 


this reason, the author elected to perform an additional 
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series of runs which incorporated a detailed experimental 
procedure. The following tasks, which compose this proce- 
dure, are discussed in their approximate order of 
accomplishment. 

The air speed and U-tube manometers were carefully set 
to zero, and all pressure line connections were checked 
for tightness. 

The Nichrome ribbon wrappings on the test cylinder 
were cleaned before each run with acetone to remove any 
film of dust, oil or oxidation. The cylinder was oriented 
into its initial position and the accesses to the test 
section were secured. 

, All electronic equipment was turned on to begin warming 
up. The initial voltmeter reading was recorded for the 
main circuit with no current flowing. This reading was 
invariably -0.065 volts. 

The wind tunnel was put in the proper gear, turned on 
and the air speed increased to only about ten mph. This 
circulated and brought the air to a uniform temperature. 

The zero was set on the millivolt potentiometer and 
the initial readings for all thermocouples were recorded 
including the initial air temperature. Just prior to 
recording these and all subsequent temperature data, a 
check was made to insure that a proper water-ice interface 
existed at the ice junction. 

The air speed was increased to provide the Reynolds 


number desired for the run. Choosing the correct air speed 
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for this purpose was not obvious. The air temperature in 
the tunnel generally increased with time to some level, 

at which the time rate of temperature rise was a minimum. 
For higher velocities the air temperature increased to 
higher levels. Since the viscosity and density of the air 
are dependent upon the air temperature, prior knowledge of 
the resultant level of air temperature was required for 
setting the air speed of the run. 

The output leads from the thermocouple rotary switch 
were removed from the potentiometer and connected to the 
chart recorder. The voltage to the main circuit was then 
increased gradually to the desired power level for the run. 
This level was chosen to be the amount of power that would 
maximize the temperature increase of the cylinder auitrace 
above the ambient air temperature without surpassing the 
“high temperature endurance limit of the acrylic tube. 
Maximizing delta T enhanced the degree of accuracy of the 
heat transfer data and tended to alleviate many of the 
surface irregularities by driving more thermal energy into 
the sealant, thereby expanding some of the recesses. 

Power was increased to each guard heater circuit until 
an approximately equivalent amount of current was flowing 
as in the main circuit; the voltage was about an order of 
magnitude less than the main circuit voltage. As the main 
eircuit settled to the desired level of power, fine adjust-— 
ments were made to the guard heaters. Their surface temper- 


ature should be maintained equal to the surface temperature 
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at the edges of the main circuit by comparing the guard 
heater thermocouple readings to those on the edges of the 
main circuit. 

After the thermal capacity of the acrylic tube and its 
foam rubber filler were saturated, a steady state condition 
existed. Then, all of the heat added was removed by convec- 
tion, less any losses due to thermal radiation to the tun- 
nel walls or circumferential conduction through the 
Nichrome ribbon. The progress of the surface temperatures 
toward their steady state values was observed and tracked 
on the chart recorder. It generally required about two 
hours to reach steady state, depending on the velocity of 
the air. The higher speed runs required up to three hours 
and the lower speed runs were sometimes steady after 1.5 
hours. 

The experiment was modeled on TRUMP, a computer program 
developed at Lawrence Radiation Laboratory, Livermore, 
Calif., for solving transient and steady-state temperature 
distributions in multidimensional systems. Reference 43 
is the manual for applying TRUMP. The time required to 
reach steady state was predicted reasonably well using the 
TRUMP model. 

After the cylinder was in a steady state condition, the 
thermocouple rotary switch output leads were again connect- 
ed to the potentiometer. The potentiometer was again 
mechanically and electrically set to zero, as it was just 


prior to any recording of temperature data, and the steady 
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state values of the surface temperature increases above 

the ambient air were recorded for the initial position of 
the cylinder. The ambient air temperature was also record- 
ed. The cylinder surface pressure was measured by the U- 
tube manometers and recorded. 

The cylinder was then rotated five degrees to a new 
position. The time required to reach steady state at the 
new position was generally 15 to 30 minutes. This was more 
than adequate time to check the ice junction, collect the 
pressure data at this new position and zero the potentio- 
meter. Once steady state was reached at the new position, 
the temperature data were recorded as before, and the 
cylinder again rotated five degrees. 

In this way, temperature and pressure data could be 
obtained at any desired angular location. All data were 
eollected on one side of the cylinder only. The experimen- 
tal runs required 8-10 hours each due to the large thermal 
capacity of the cylinder, and the small number of test 


thermocouples installed. 
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VI. PRESENTATION OF RESULTS 


Raw data, including temperature information in milli- 
volts and pressure information in cm or in of water, were 
read directly into the computer program in Appendix B. 
The density, thermal conductivity and kinematic viscosity 
of air were evaluated at the film temperature, Tp» which 


was computed by 


(T (ae 


alereia a 


ani! avs f 0.5(AT) ove (43) 


) is the average ambient tem- 


for each run, where Cres ane 


perature and Ee is the average local surface tempera- 


ture increase at steady state. 


A. TABULAR RESULTS 

The computer program provided detailed results for all 
thirteen experimental runs. These are presented in tabular 
form on pages 72-84. The pressure coefficient, heat trans- 
fer coefficient, Nusselt number and Froessling number are 
provided at five deg intervals from the forward to the rear 
stagnation points for all runs except three runs which 
omitted the pressure information. Two of these runs are 
duplicates of others which include pressure information. 
Therefore, a complete set of results is presented at ten 
different Reynolds numbers. 

"DELTAX" is the are length, "ASURF" is the surface area, 


and "QGEN" is the heat generation rate of a 5 deg element 
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of Nichrome ribbon. "TU" is the turbulence intensity level 
at the given free stream velocity. "DELTA-T CORRECTION" 

is the correction for the resultant effect of viscous heat— 
ing of the cylinder and heating of the ambient air during 


high speed runs (Appendix C). 


B. DISCUSSION 

To illustrate the discussion, the Nusselt number and 
the pressure coefficient have been plotted as a function of 
their angular location in Figures 9 through 14. 

The results in general only confirmed the trend of 
those published by Giedt [15] with regard to approximate 
locations of minimum and maximum points and the angular 
pattern of heat transfer characteristics around the cylin- 
der. The magnitudes of the heat transfer rates are some- 
what less than those of Giedt. This is most probably 
explained by the lower turbulence level of the wind tunnel 
employed in this experiment (~0.5%). Also, the alignment 
of the cylinder was about one deg in error, as it leaned 
toward the rear of the test section. This discrepancy 
probably decreased the normal component of the wind strik- 
ape ne cylinder, and may have ‘consistently decreased the 
heat transfer rates. 

At the lower Reynoids numbers of 57,000 and 82,000, the 
heat transfer results were in close agreement with the 
resulting pressure coefficient distribution around the 


cylinder. The boundary layer remained laminar and 
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Figure 9. Heat Transfer Results at Reynolds Numbers of 
57,000, 223,000 and 350,000 
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Figure 10. Heat Transfer Results at Reynolds Numbers of 
82,000, 261,000 and 406,000 


87 


co) ae ‘ 
ee ie . er an _ 









000,$6 = sa 0 
t= 060,185 = of 2 | 
000,308 = 474 










sa ey Ot! t —_ 2 _ ~ es | Als b= my 
31, Gd eB. Oa OOL os Gg ov og G 


sisdmu4% eGfonys t e2?fue ~plenen? gs0H ,0L © 
» bre 000,108 (000,58 


- 
a 






700 


600 


Nusselt Nr... 


300 


200 


100 


Figure 1l. 





o Re = 153,000 
0 Re = 309,000 
= 495,000 
LA 
is AS AA a 
ALLA 
6,6,8 
2 7 o 
oO 
WM—-5 
oO OOO O 
Q 
% 1) 
OL, O. 0 oe 
O Oo 


20 4O 60 80 LOO: 120 AAO’ ~ LEQ ~ USO 
Theta 


Heat Transfer Results at Reynolds Numbers of 
153,000, 309,000 and 495,000 
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separated at about 82.5 and 80 deg, respectively. This 
laminar flow regime was reflected in the heat transfer 
results by a single minimum value of the Nusselt number, 
which was located at the separation point. Also, maximum 
heat transfer occurred at the stagnation points. The cor- 
responding pressure distributions also depicted laminar 
flow. These were characterized by a minimum at 65 deg 
followed by a fairly small recovery in the form of a 
smooth, continuous rise, and a flattening out at the sepa- 
ration points. Additionally, there were distinct peaks in 
the two pressure coefficient curves at 85 and 80 deg, res-. 
pectively, just as the curve flattened out. 

Unfortunately, at the intermediate range of Reynolds 
numbers studied (82,000<Re<300,000), the surface pressure 
distribution did not reflect the same conditions of flow 
inferred by the heat transfer results. As before, the 
pressure coefficient distributions were representative of 
subcritical flow, and indicated a laminar boundary layer 
separation at about 80-85 deg. The heat transfer results, 
howéver, depicted a mention One cranswie tonal ivowlsn tne 
boundary layer, and therefore were representative of crit-— 
ical flow. There were, for example, two minimum values of 
Nusselt number. The first minimum occurred at an angular 
location of about 85 deg, and marked the beginning of 
transition to a turbulent boundary layer. At this point 
the boundary layer may have separated, become turbulent, 


and then reattached at a subsequent point, having become 
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fully turbulent. This transition region, known as the 
"Separation bubble", is exhibited by the angular distribu- 
tion of the local Nusselt number, at this range of Reynolds 
numbers. The second heat transfer minimum occurred at about 
130 deg, and represented the separation of the turbulent 
boundary layer. At this range of Reynolds numbers, and 
between the two minima described, maximum heat transfer 
occurred at an angular position corresponding to the re- 
attachment of a fully turbulent boundary layer. 

The lack of agreement of the pressure information can 
most probably be attributed to the location of the pressure 
pickups on the smooth acrylic tube and not in the region of 
the test section. 

At the higher Reynolds numbers of 406,000 and 495,000 
there was close agreement between the heat transfer and 
surface pressure results. Both indicated critical flow 
conditions. For each run the heat transfer results exhib- 
ited two minimum values of the Nusselt number. The first 
minimum values where transitional flow began were located 
at 90 and 82.5 deg, respectively. The second minimum 
values where the boundary layer separated were located at 
130 and 135 deg, respectively for the two runs. The points 
of maximum heat transfer where the boundary layers reat-— 
tached and became fully turbulent were located at 100 and 
95 deg, respectively. 

For these same two Reynolds numbers the surface pres- 


sure distributions had minimum values at 80 deg, which 
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differed from all other experimental runs which had pressure 
minimums at 65 deg. Also, the rise of the pressure curve, 
from the minimum, was not smooth, but was characterized 
by a double inflection, which depicted the "separation 
bubble." The first inflection point represented the initial 
detachment of the laminar boundary layer. Its angular loca- 
tions of 90 deg and 85 deg for the two runs corresponded 
identically with the locations of the first minimum values 
of the Nusselt number. The second inflection point indi- 
cates reattachment and complete transition to a turbulent 
- poundary LACE eriosSmaneilar Locabionmon l00F dep am: both 
runs, corresponded closely with the location of the maximum 
values of the Nusselt number. Also the pressure distribu- 
tions for these Reynolds numbers flattened out at angular 
locations of about 130 and 120 deg, respectively, corres-— 
ponding reasonably well with the second minimum values of 
the Nusselt number. 

In contrast to Giedt's findings [15] for Reynolds 


numbers above 153,000, the value of h, at the rear stagna— 


8 
tion point was not always greater than its value at the 
forward stagnation IOOTING 

It was not found that the angular location of minimum 


values of h, invariably increased as the Reynolds number 


8 
increased. Instead, there were variations in the flow con- 
ditions from one run to the next. A higher velocity run 

at times produced an apparently smoother flow than one of 


the lower velocity runs. For example, the initial minimum 
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values of Ng for Reynolds numbers of 153,000; 173,000; and 
223,000 were 90,100 and 100 deg, respectively. However, 
at a Reynolds number of 260,000, the first heat transfer 
minimum occurred at 82.5 deg. In the same way, the initial 


minimum value of h, for Reynolds numbers of 309,000; 


8 
350,000 and 406,000 were 85, 95 and 90 deg respectively. 
But apparently a very smooth flow existed at a Reynolds 
number of 495,000, since the first heat transfer minimum 
then occurred at 85 deg. This inconsistancy was probably 
caused by the surface defects of the test cylinder. These 
defects were less pronounced during the two high speed 
runs because of the higher surface temperatures involved. 
It is thought that the large thermal capacity of the 
foam rubber and acrylic in the interior of the tube fos- 
tered complicated internal flows of heat. These heat flows 
depended upon the magnitude and direction of the angular 
rotation of the cylinder from a previous steady state posi- 
GOMmGO Uhe present position. Also, the magnitudes of the 
steady state surface temperature differences varied, depend- 
ing upon whether the cylinder was rotated to a position 
from a previous steady state position, or was initially 
warmed up there and heated to a steady state condition. 
Specifically, and for example, rotating a test thermo- 
couple at five deg intervals gradually around a distance 
of 20 deg to a critical location seemed to amplify these 
temperature differences than if the cylinder had initially 


been warmed up with a test thermocouple in that critical 
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location. This was especially true when approaching and 
arriving at the first minimum value of the Nusselt number. 
Variations near the front of the cylinder were only about 
0.5°F or less, but were on the order of five °F at the 
critical position, and about one °F at the rear stagnation 
point. 

An uncertainty analysis was performed and is described 
in Appendix F. In this analysis the uncertainty of the 
Nusselt number was found to be about three per cent at the 
forward stagnation point, about six per cent at the rear ~ 
stagnation point, and on the order of 25% at the critical 
point where the boundary layer either begins its transition 


to turbulent flow or separates. 


C. COMPARISON OF EXPERIMENTAL RESULTS WITH EXISTING THEORY 

The heat transfer at the forward stagnation point can 
be conveniently compared to Squire's solution, which gener- 
ates a Froessling number (F = Nu/YRe) of 0.995 for this 
application. The Froessling number calculated from the 
present set of experimental results was about 0.93-.92 at 
Reynolds numbers from 153,000 to 261,000. At Reynolds 
numbers below and above this range, the Froessling number 
decreased to as low as .84, excluding the run at a Reynolds 
number of 495,000. This latter run had a surprising, low 
Froessling number of 0.76 at the forward stagnation point, 
which could not be explained. For this run the power 


source for the main circuit and the wind tunnel were both 
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in their maximum positions, which possibly placed an unac- 
ceptable strain on the experiment. It is also curious to 
note that the Nusselt number at the forward and rear stag- 
nation points are identical for the two experimental runs 
of 406,000 and 495,000. 

For an evaluation of the data beyond the forward stag- 
nation point an appropriate theoretical solution was re- 
quired. All existing theoretical solutions only predict 
the heat transfer from the forward half of a cylinder, 
prior to the point of laminar boundary layer separation at 
about 80 deg. Spalding and Pun [12] have summarized free 
teen approximate solution methods. However, since all fif- 
teen are applicable to isothermal cylinders, none would 
provade a meaningful comparison for the data collected in 
this experiment with a constant heat flux cylinder. 
Spalding and Pun compared the various methods with that of 
Froessling who used an "exact", series method of solution 
limited by the number of terms employed. It was shown that 
Schuh's method of solution for an isothermal eylinder 
ranked in the highest class of accuracy and within 1-3 per 
eent of Froessling's solution. Schuh also published [26] 
one of the few approximate methods of solution for a con- 
stant heat flux cylinder. 

. In Fig. 15 the two theoretical solutions of Schuh are 
DleGtced to filustrave the significant difference in the 
heat transfer profile between the two types of cylinders 


usually employed. As one would expect, the two solutions 
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converge at the forward stagnation point, for in that vicin- 
ity the constant heat flux cylinder also has a constant wall 
temperature. Moving farther around the cylinder the ratio 
of relatively greater heat transfer from the constant heat 
flux cylinder continually increases. This is physically 
represented by an increasing surface temperature, which is 
required to provide an increasing temperature difference 
across the thermal boundary layer. The increased tempera- 
ture difference is required to drive a constant flux of 
thermal energy through the thickening thermal boundary 
layer. 

Some of the experimental results of Giedt [15] for 
heat transfer from a constant heat flux cylinder are com- 
pared to Schuh's solution in Fig. 16. The nature of the 
divergence of Giedt's results from the theoretical curve 
is most probably explained by the significantly higher level 
of free stream turbulence (2.25% compared to 0.5%). 
In Fig. 17 some of the present results are compared 
with Schuh's solution, and most of them agree within about 
six per cent. This agreement is largely due to the low 
turbulence level of the wind tunnel used (0.5%). These 
results indicate that slightly less heat transfer occurs 
than that predicted by Schuh. The results are closely 
grouped, implying that any experimental errors are fairly 


consistent. 
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Figure 16. Comparison of Giedt's Experimental Results to 
the Theoretical Solution of Schuh 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


Heat transfer and pressure coefficients were determined 
around a heated, right circular cylinder placed in a cross 
flow of air. The heat transfer results were generally in 
agreement with the pattern of heat transfer found previ- 
ously for subcritical and critical flow. 

The heat transfer results were compared with the clas- 
sic results of Giedt [15], and indicated relatively lower 
heat transfer rates. This is most probably explained by 
the considerably lower free stream turbulence level of the 
wind tunnel employed in this experiment. It has been pre- 
viously shown by other researchers that higher free stream 
turbulence levels markedly increase the heat transfer 
rates. 

The heat transfer results were also compared and showed 
good agreement with the approximate analytical solution of 
Schuh for the laminar flow region om Chee preonita porto Ot 
the cylinder. 

The surface pressure distribution did not in all ex- 
perimental runs support the heat transfer results. Specif-— 
tieally, in several runs the pressure information failed to 
indicate the "separation bubble" which represents the trans- 
ition in the boundary layer to turbulent flow, and depicted 
subcritical flow instead. In these instances the heat 


transfer results clearly represented critical flow effects. 


102 





nal! - 


ined  eaiead ylovidatet betanthnat brs. ,[2L] ta 0 to atts 



















bentereseb 91184 ‘sinstoftteoo oteeena Bae setanexs § 
380% a ai beoata rabatiys zekuaite’ “ata ky cbedead, a 
pa el tircenas wrsW esluaet astenses taem ‘aa ak 
 ntvetq brvot ws'tenpsd+taet to aeseq, ote ‘dgkw ta 
swolt ISOS EED dne-inold trode 101: 


Td aaa ie we er 


<aals ord cdolw DHS UNOS ore ‘edtiuses is inne es : 


vd bentsliqzea yidedorwq Jaom el slat A987 s9tse 
sdPIS’ isevel sonaludw2 mserte sat tewol recall 
-8tq ased ead JI pee efad ad bayolqms te 
mesise av inaety be tart? a4 edowscest 1Sdio xd words) 

selene ty Fasc sy sessteat <ibedasm elevat oF 
| “ef 


bewotls boa betsqmon osie etow ediuest telaneyd Jeod em 


10 noliuloa isolsulens edamixotags odd déiw: daeneata 


aA, 


Ye soloed Jnott edt no oolges wolt tanigal aad cee 


tig 


ke Iie ut Jog BIb notsodiagsthth ST ReoTq | aoe les | 


~lioeqe .edluesey 79ers 13 Jead ef Itogqua amis ising vy 
oo heltsel solgsmiotnl siwanetq of3 save Leteves ah ue tat 
<“niien? ons evieaptqer told “sled netieageee was o2u0 Bt [ 
befolgeb bra welt Jsnoluduus of sogal yrebawed oi2 ab nods 

afd sienssent ‘Seodd aI J bsetant wot} CAo13 trodua. 


8150TTo wold Capt? tio Davnsesizqes ¢{itsels wiiveor tetanes 


sob’ & man 6 


It is possible that the pressure measuring system was not 
adequately sensitive to detect all of the surface pressure 
variations. However, it is more probable that the loca- 
tion of the surface pressure pickup was too distant from 
the test section. It is recommended that on the next cyl- 
inder constructed, the pressure taps be located in the test 
region. 

An uncertainty analysis was performed and showed that 
the uncertainty in the heat transfer results was about 
three per cent at the forward stagnation point, six per 
cent at the rear stagnation point and on the order of 25% 
at the peaks in the heat transfer curves, where fluctuations 
seemed to be amplified. 

The heat transfer results were possible adversely 
affected by some defects in the test cylinder, which can 
be avoided in a second cylinder being constructed. 

The easily deformable acrylic tube and its large ther- 
mal expansion coefficient as compared to Nichrome ribbon, 
and the curing properties of the sealant used to bond the 
Niehrome to the surface resulted in a number of surface 
irregularities. A plan for avoiding this defect is being 
implemented in the construction of the second cylinder. In 
the "Discussion" by B.H. Sage appended to Ref. 32, small 
surface irregularities are reported to have caused unusual- 
ly large variations in the heat transfer, especially in 


the rear portion of a cylinder. Seban [32] stated that 
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similar surface defects "disturbed the flow near the separa- 
eatoyg Tstigves 1 

It is recommended that the TRUMP computer program [43] 
be utilized to evaluate various insulating materials as a 
suitable replacement for the foam rubber. It should be 
possible to eliminate the complicated thermal flows and 
rotational behaviour, long transient warm-up period, and 
possible internal free convection as the acrylic tube 
usually expanded clear of the foam rubber filler. 

It is recommended that the cylinder vertical alignment 
be carefully checked before all experimental runs. The 
base plate moving parts and metal interfaces should be 
overhauled, sanded smooth and lubricated. The rotating 
dise should be marked to indicate every one deg; it is 
the author's opinion that in a tunnel with a low turbulence 
levels a sufficiently accurate study is possible that 
Wakroancs this rerinement. The wooden floor has become 
warped from heat and should be renewed. 

A major source of concern with regard to accurately 
measuring the Bree Gt pean velocity was the air speed mano- 
meter. It is recommended that a hot wire anemometer be 
installed near the cylinder in such a manner as not to 
disturb the flow. This also precludes the problem of in- 
volving corrections for wind tunnel boundary constraints. 

It is also recommended that a separate study of the 
same cylinder be made using cholesteric liquid crystals as 


a temperature sensor on the surface of the cylinder. 


104 


sessase wit sai wort) ond easter 


PS as a 


tee A 
Oe | ‘any ; 
nh iy = Ory oe oe i 


ri eee oe ane _ ' | oe At ‘ ; j 
iN .) LY Db eee 5 I . me : Uo 7. 2 baie ‘ i ' ve 7 
Rael Tt Be Crny . oy ae a 






ra) ey Ste un i Ps iz 

: ; thy i ie one oe 4 » Bs ba A, SR : aien , £ one : 
“tei, margotg oduqmos aucun ae Sata; Feb: abe ahs, 
ps as elehrotein aitistvent amobiay eduglevi’ 0 beneta 7 


od bivora, ¢T , -naddues meet ane 202 oe 
ban awolt Leora? batsoliqnys _ otentativ ot eb 
























bas ,botteq qu-m iw gosténatt anol ,wolvadsd sine 
edud oflytos gid ea Holtpevaog oef?-Lanistab sidfi 
«I0 L527 4oddua maot odd ds tes Hahoaqne 
sammageta feoluusy tebdh tvs o.l9 Jang bebabmnose4 att 
sdf .znuy Lsinemtisqro ifs stoted. betsedo yl 
od bluors-eenatiooni Letvem baa eftaqd yalvom Pey 4 
aAatisiot oe «beJaolidul bas dsoome bohase-bsluad 
el @t 4348 “to yueve e¢solbnt b4°héediem od bf idl 
eonéeiudiuy wol 2 djiw fennud s ni gané oolaioqa a" totsis 
Jai? sidlenog a. ybwie etsuuoss ylinatolht ives ie 
omosed 26d tool) meboow ed? .duementies atte 2 
«fHewonet- sd Sluods Sha seed modi pst 
Yistpusss oo brage: ddiw: miooreo to santos cotems 
-onam 5Bsege tis sid saw yttoclsv meetia ast? edd wclige 
od aeJomonsns. o1tw tod ¢ Jane bshashmoset al oT ne 
ov VO @s Jenirss s dogs al sebantive BSA? ts98 Hel val 
at Yo mofdotq ens sebuloe xq oslé sldT .wel? on? dwitel 
BILE Bro® ershMod Jeans balw io anolase'rmes ant Low 
oH4elo “ebudeedaieges a tad? bedriommoned Waar aL lee 
BB fESINY IO phagy 5 aossateda antow eham ad Setriiys smag | 
“ebative i? to sostigs oft AS soenesS st Faro qed s 


- 


These crystals indicate temperature through a change in 
color, and offer the possibility of providing a highly 
desirable visual display of the temperature field much as 
smoke streaks provide a visual display of various flow 


phenomena. 
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APPENDIX A 


Wind Tunnel Blockage Corrections 


The test section velocity obtained from the air speed 
manometer and its associated pressure pickups refers to a 
clear test section. In Ref. 42, Pope discusses the effects 
of the tunnel boundary restrictions in combination with a 
two-dimensional model mounted in the test section. MThe 
flow area is reduced by the physical presence of the model 
and its relatively slower moving wake. Both of these dis-' 
turbances cause the air to accelerate as it flows past the 
model in order to preserve continuity. Hence the measured 
velocity must be corrected for solid blocking and wake 


blocking, respectively. 


Calculation of the solid blockage correction: 


This boundary correction was calculated using the 


equation of Pope Ca 


Coie) & ‘es gt (44) 
en) 5) | 

on 4 [Ge] 

Egan = 0-008215 


Calculation of the wake blockage correction: 


This boundary correction is a function of the drag 


coefficient, and had to be evaluated for each run. The 
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equation used is from Pope [42], and the drag coefficients were 


were taken from Fig. 1.4 of Ref. 5. 


WB (+ (2) co (45) 
“WB” = 
Run ene es ee WB te OSB) CWE 
1 54,785 1.20 0.03 0.038 
2 79,296 120 0.03 0.038 
3 97,522 1-20 0.03 0.038 
y 147,675 120 0.03 0.038 
5 166 , 381 1,520 0.03 0.038 
6 215,044 120 0.03 0.038 
7 215,185 1.20 0.03 0.038 
8 250,655 1.05 0.026 0.034 
9 251,928 1.05 0.026 0.034 
10 298,958 0.95 0.023 0.031 
wa 340,095: 0.85 0.021 0.029 
ne 305,087 — Ono Oo Oly 0.025 


13 487,288 0.30 0.007 OROM5 


Application of the total blockage correction 
The right column of the table provides the total block- 


age correction required for this experiment. Only the free 
stream velocity and the Reynolds number were corrected for 


boundary effects. The boundary corrections did not enter 
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into calculation of the pressure coefficients. 
Specifically, 
U corrected = U uncorrected (1+e,) 


Re corrected = Re uncorrected (1te,) 
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Computer Program for Data Reduction 
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APPENDIX C 
Wind Tunnel Calibration for the Combined Effect 
of Ambient Air Temperature Rise and Viscous Heating 


of the Cylinder 


In general the free stream air temperature steadily 
increased during the operation of the wind tunnel. This 
temperature rise was attributed to the viscous heating in 
the boundary layer of the wind tunnel walls and the energy 
added by the rotating impeller. At low speeds (Re<260,000) 
this temperature rise was very gradual and did not affect 
the accurate measurement of AT. The measured AT represented 
the temperature rise solely due to the applied power. Al- 
though the surface of the cylinder became warmer due to the 
heabed yair this extraneous heating of both the cylinder and 
Ghescleewere at the same rate, and AT remained a true 
measurement. During high speed runs, however, the rate 
of temperature rise of the tunnel air was rapid and ex- 
ceeded the corresponding cylinder surface temperature rise. 
Ths effect caused the measured AT to be smaller than it 
should have been. 

AG the same time, and also for the high speed runs, the 
viscous heating of the cylinder surface provided a counter 
epmheet, tending to increase the measured AT. In all of 
the high speed runs this latter effect was dominated by the 


ambient air temperature rise. 


ALL) 











O RIGUSIGA A 
) 4 
doors bontamss oat od aobinwétad fs 


enn 


yatiasy auooetv Mid eult’ uy Lone ao ak anol PO 


= os: 


. , 
; ‘i q 
* ton ® cunern sd¥ to " ys 


em» MEtbaode olidertageset Kee une sert sd? Coons 
| bidy. sf annnd. barby ori 3c. naltqrege. ‘ote gntnwd ba 
“pb anddestt auogely att ba bedudtntis ne, sain 

| | Yates sit bie sifow Loniurd btw ord 4a “oyel 
(O00, 005298) Sbesqs wol JA isnifegnt antsason ae 


Joetta Jon 216 Oris, levbara yvxev sew oelt 









| aa 
an 


_betecsorget Th benvetem sxT. .TA 6 Taomsrnlangt Ye 
fA 18Wed bolfaus edt of seb) vlelo&® seks Deke 
93°69 sub. roui'taW emsoed tshaltlys ong to spats ony! 
bis tebniive sny.diod to griveed sucshanese etal 
outt 2 hentamoy TA Sos , stat mek odd 3s 7s 


2S S02 ,tevewon , 4 beags dala gntagd 79 











_— 


~xe Ota Oiqes sew tla fennud ets Yo Sekt Souves 
9829 sioJnteqmes sob tive teaballye sotbnogas?i0s saz 


$i ted? tolfane od of TA botvbaen 449 Deeuar JoaTls 

| Wrend Svea | 

sid , entry Sesge dali shY tot coffe bra oak ofiss. sag shaq 
tS7nNvOo B Seblvom sostavs yeintiive out to nuivsed eavos 

2 | 20 digs al ..TA Dextesgemiedd me aotant ‘og anliies 408 : 
MSE YS Sotectomh agv soette totjal side aus hesqe gta onty 


-onln or tena qed ate gastome 


Sit 


Rather than attempting to uncouple these two effects, 
a scheme was devised to calibrate their combined effect 
CULT On Nunc dat ay forthe rhtehwispeedi runs | lhe) Natchrome aati — 
bon was not energized and the cylinder was heated only by 
the extraneous effects discussed in the preceding para- 
graphs. After several hours steady state conditions were 
determined using the same technique as described in Section 
V. EXPERIMENTAL PROCEDURE for an ordinary run. The temper- 
auuGe ian herence (Tp-T oar) was then measured at various 
angular locations around the cylinder, and treated as a 
correction factor for the data collected during an actual 
run with applied power at that Reynolds number. 

Calibration runs were performed at Reynolds numbers of 
309,000, 406,000, and 495,000, and the AT correction factors 
are provided in Section VI-A TABULAR RESULTS. In the data 
reduction computer program (Appendix B) the AT correction 
factors are subtracted from the measured AT just prior to 
the calculation of the local heat transfer coefficient. 

The radiation and conduction losses are based upon the un- 
corrected, measured AT. The AT correction factors are most 


often negative because the heated air effect was dominant. 
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APPENDIX D 


Wind Tunnel Turbulence Level Measurement 


The turbulence level of the Acrolab wind tunnel test 
section was measured with a Thermo-Systems hot wire anemo- 
meter at eleven different free stream velocities ranging 
from 10 to 110 mph. 

Initially, a Thermo-Systems calibrator model 1125 was 
used to calibrate the probe. The cold wire (room tempera- 
ture) resistance was 7.77 ohms. With zero air flow through 
the calibrator the initial probe voltage was 1.667 volts. 
ihewprobe voltage was then recorded for various controlled 
values or air flow velocity through the calibrator. The 
results of this calibration are presented in Fig. 18, as 
aeplot of the probe voltage squared, View versus the square 
root of the air flow velocity through the calibrator, VU. 
ine intercept of the linear calibration curve anaien result- 
ed, was the square of the actual value of initial probe 


voltage, V to be used in the turbulence level 


92 
calculations. 
The probe was then mounted in the center of the test 

section on a vertical rod whose height could be adjusted 
in one inch increments. At each air speed.for which turbu- 
lence level measurements were taken, the rod was used to 


traverse the vertical height of the test section. The 


probe voltage, V, and the rms voltage, (ec, were recorded at 
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Note: V2? = We + BYU 


Calibration Curve for the Hot Wire 
Anemometer 
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were observed to be constant for any given air speed except 
near the floor and ceiling where wall boundary layer ef- 
fects caused some fluctuation. 

Themouuchuwlence level Lu was Calelledteduca tut Mem mjc OUs 


air speeds by making use of the relation 








where U represents the mean free stream velocity, and /y'? is 
the mean fluctuation in the free stream velocity. The 
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APPENDIX E 


Endurance Test for the Acrylic Tube 


An eighteen in section of the 4.5 in o.d. cast acrylic 
tubing was wrapped with six turns of the Nichrome ribbon. 
The ribbon was instrumented with one copper-constantan 
thermocouple. The surface temperature was then measured 
using a Leeds and Northrup millivolt potentiometer and an 
ice bath for a reference. A Lambda 60 volt power supply 
was used to provide a 2.0-2.5 amp current through the 
ribbon. The tube was thereby heated, uniformly and slowly, 
TO CeMmperauune shevetismot OOH. LO. a 20° vand 50° fh, 
and cooled slowly re room temperature after each heating 
cycle. The heating and cooling rates were only about one 
to two degrees per min., to minimize the relative differ— 
ences between the acrylic and the Nichrome during the ex- 
pansion and contraction. In all of the heating cycles 
excepu the final one the cylinder did nov physically de- 
form, as there was not a noticeable expansion of the 
acrylic. Tm the final test, between 130-135°F the acrylic 
began to visibly expand causing the ribbon to tighten. 
This immediately caused a hoop stress on the surface which 
deformed the cylinder by constraining its expansion in the 
area where the ribbon was wrapped. The cylinder expanded 


further in the unwrapped areas. When the cylinder cooled 
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it contracted uniformly leaving a smaller, deformed diameter 
ineche recon lon. the wlbbien. 

Rrom this test it was concilluded that to vavoid the risk 
of deformation the actual test cylinder should not be 
heated above 120°F. 

In a second test the surface temperature was heated to 
180°F, severely deforming the cylinder. However, the 


R.T.V. sealant was unaffected and held solidly. 
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APPENDIX F 


Uncertainty Analysis 


ASPan cillnexperimMenuak SivUdLes.. Mere existed ssOme 
degree of uncertainty with regard to the accuracy of vari- 
ous results, which stemmed from the data collected. To 
evaluate the levels of uncertainty involved, the method of 
Kline and McClintock described in Ref. 44 was employed. 

The measured variables which were the origins of the 


AP AT, V, and R (see no- 


uncertainty were INE ; cylinder? 


menclature for the definitions of these terms). The fluid 
properties of air were extracted from table A-7 of Ref. 45, 
and their uncertainty is neglected in the analysis. 

The magnitude of the uncertainty associated with the 
pressure measurements depended upon the flow velocity of 
the run, and the angular location on the cylinder. Simi- 
larly, the magnitude of the uncertainty associated with the 
temperature measurements depended upon the angular location. 

In this appendix the calculations of uncertainty are 
Snowe in detail for a representative Reynolds number of 
260,600, and at the forward stagnation point of the cylin- 
_ der. Uncertainties for other flow velocities and locations 


are discussed. 


The Pressure Coefficient 
The pressure coefficient was calculated from equation 


28, 
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Applying the technique referenced above’, one has 


2 W 2 
Bender ‘ APi—2 
AP,_-» g 


| cylinder 





where WwW is the uncertainty interval, or precision range. 


At the forward stagnation point, 


= + 
AF oylinder (Goa On eMmewacerGUO ston 4) 


een Onlencm: wawer Gl. to) a) 


ME 


Substituting, one has 





inerenore nunc UNCertalnby in the pressure coefficient was 
about two per cent. In the wake, the uncertainty was about 
four per cent because the unsteady motion of the air les- 
sened the readability of the manometer. 

At the higher Reynolds numbers the uncertainty in the 
pressure coefficient was about three per cent at the for- 
nae SGaenehaen point and about six per cent in the wake. 
The higher speed flows produced larger fluctuations on the 
manometer. 

At the lower Reynolds numbers the measured pressures 
were small and the uncertainty intervals were a significant 


fraction of these pressures. As a result the uncertainty 
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in the pressure coefficient was large, reaching as high as 


49% for the lowest run. (Re = 57,000) 


The Heat Transfer Coefficient 
The heat transfer coefficient was calculated from the 


equation, 








where A. is a constant (see nomenclature for definition of 


terms). Continuing, 


W QW 2 wW ; W ‘ 
J aL eas + iss + _AT 
h Vv R AT 


For the experimental run at a Reynolds number of 260,600, 


FOO ae nv OMeticn Gl Ontion a) 


< 
i] 


mw 
i) 


20a One molm si @h Ont. ones) 
ME A000) = On5ch (lO to 1), at the forward stagnation 


point. 


Substituting, one has, 


Ww 2 2 2 

ma A s OFS 0.2 ORS 2 

h [2 aa i [38-2] : [so:3| me Gus : 
The uncertainty in the heat transfer coefficient was about 
three per cent at the forward stagnation point for all the 
experimental runs. 


At the transition point where the boundary layer be- 


comes turbulent, or at the separation point in laminar 
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flow, the temperature difference had an uncertainty inter- 
Valcom abou veon and swhesuncertiainty in the heat trans— 
fer coefficient was on the order of 25%. 

The uncertainty of the heat transfer coefficient in 


the wake was about six per cent. 


The Nusselt. Number 


The Nusselt number was calculated from equation 3, 


Continuing, 


and 


i) NN se OS) alin (ClOmtsion al)ys 


The wuncervainty of the measurement of the diameter was 
ne wMeaENCh mE inere hOre une Uncertainty im the Nusselt 
number was the same as the uncertainty in the heat transfer 


coefficient. 


The Reynolds Number 


The Reynolds number was calculated from equation l, 


and its uncertainty is the same as that of the free stream 


velocity, U, which was calculated from equation 29, 


ad Ss 
p 


723 


cm) ad 7 
RPO aa A a OR 


onsen rw 

















‘2s 3 ie earnoiadt Parcies ion 


: st: ot dndtolte02, “orenstt taedt add to mee son 

—— Lf a 
Da hee eae vps 709 ata eds iw ol 
i‘ : 


. | ae, es 
nolspupe mow? Desefwolas sew nsdmuh diseer 


as | 





abe fe @ : | 


-(L'od 32) ah 20, + ye. 


\ Bow sstemaih af7 To snemeiwepet Sid 30 Yonteeteee 
aS 
tloeag ef ol vinistresnd sad ,saaieyen? ~ sige 
ow 

: ~ well 
7 toling1s Jaen sid nl ysyAtayssody SAF GA SmAS end Sane 
~ ~ 

ain a 

7  Inets. 


, 
r 


rs 7 e 
SedaeK sh Lomge 
; i al 


.f nofzaips wout besetustsa-aaw togmun sbionyenta 


riaitea oeTl gif 14 Jany 66-emar acd ak SONA e 50H oat Bi 





; eS nol ’ tatueteo Baw dotely yy <a 008g 
>. (gas 
; q 


Continuing, 


WW) W 3 : 
RES creUis a s 
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Therefore, the uncertainty in the Reynolds number is approx- 
imately half the uncertainty in the dynamic pressure, and 
equals about 0.7%. At the high Reynolds numbers this un- 


certainty doubles. 


The Froessling Number 


The Froessling number was calculated from equation 4, 


me Nua 
¥Re 
Continuing, 
\2 
OE = Me Nu WRe 
F Nu Re 





Substituting from previous results, one has 


W a 
= = (0.03)? + (0.0035)? = 0.03 

The uncertainty of the Froessling number was generally 
about the same as that of the Nusselt number. At the for- 
ward stagnation point the uncertainty of the Froessling 


number was about 3-4% in most cases. 
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